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ABSTRACT 


A hyperbaric chamber environment is simulated using a 
FORTRAN computer program. This program allows the user to 
input data files of varying chamber conditions and provides a 
detailed output which can be used for design or chamber 
operations. Specifically, this program determines life 
Support, temperature and humidity control, and power 
requirements. After developing the theory and describing the 
program in detail, complete examples simulating a hyperbaric 
chamber at 2000 fsw, a submarine lockout chamber at 400 fsw, 


and a shallow underwater habitat at 50 fsw are demonstrated. 
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CHAPTER I 


INTRODUCTION 


A hyperbaric environment exists ina facility which may 
be pressurized from one to several hundred atmospheres. 
Edmonds (1976) defines a hyperbaric chamber asS_~— an 
environmental pressure vessel where the ambient pressure can 
be increased from one atmosphere to many times atmospheric 
pressure depending on the structural capabilities of the 
vessel concerned. Hyperbaric chambers are also known as 
recompression chambers and decompression chambers. A 
recompression chamber is used for compressing a patient or 
diver as part of a therapeutic regime or for training. A 
decompression chamber is used to decompress a subject already 
exposed to increased pressure. Miller and Koblick (1984) 
define an underwater habitat as a facility which acts as a 
survival shelter allowing humans to survive and live on the 
ocean floor continuously for several weeks. Hyperbaric 
chambers, underwater habitats, diving bells, and submarine 
lock-out chambers require unique design considerations to 
account for the additional pressures and the effect on human 
physiology and thermodynamic properties. 

History 

Hyperbaric chambers have had medical uses for over 300 
years. Edmonds (1976), Miller and Koblick (1984), Shilling et 
al. (1976), and U. S. Navy Diving Gas Manual (1991) provide a 


chronological history of hyperbaric chambers: 





1662 - The first use of a hyperbaric chamber in medicine was 
made by a British physician, Dr. Henshaw. This chamber was 
pressurized and depressurized by a pair of organ bellows. 

1830's - Three Frenchmen, Junod (1834), Tabarie (1838), and 
Pravaz (1837) renewed interest in hyperbaric chambers. 
Junod's chamber was a 1.6m diameter sphere capable of 
pressures up to 4 atm. He noted the euphoric effect of 
compressed air. Tabarie claimed improvement in the treatment 
of lung disease. Pravaz built a large Ben her capable of 
holding 12 patients at a time. 

1850's - A rapid rise in the construction of hyperbaric 
chambers was observed. Compressed air was provided by steam 
engines and chambers capable of simultaneously treating up to 
50 patients were built. 

1860 - The first hyperbaric chamber in North America was 
built in Ashawa, Canada. 

1891 - The first chamber in the United States was built in 
Rochester, NY. 

1921 - Orval J. Cunningham of Kansas City, MO constructed 
the world's largest hyperbaric chamber. It was 5 stories 
high, had a diameter of 20 meters, and was equipped with 


amenities suitable to a first class hotel. 


Since 1920's - Numerous hyperbaric chambers have been 
designed for treatment of divers. With the progress of 
hyperbaric and diving medicine, chambers are becoming 


increasingly complex. 


1960's - The French and Americans began work on underwater 





habitats. The French built Conshelf I & II while the United 
States built Man-in-the-Sea I & II and Sealab. 
Modern Hyperbaric Facilities include: 

1) Ocean Simulation Facility (OSF) - located at the U. S. 
Navy Experimental Diving Unit in Panama City, FL. This 
chamber is one of the world's largest man-rated hyperbaric 
facilities with a depth capability of 2250 fsw. A wide range 
of experiments in biomedical studies, human factors 
verification, and deep diving systems performance tests are 
conducted in wet and dry chambers. 

2) Naval Medical Research Institute (NMRI) - located in 
Bethesda, MD. This complex is the U. S. Navy's largest 
hyperbaric complex fully devoted to research and contains 
three dry chambers and a vertical wetpot rated to 2280 fsw and 
two dry chambers rated to 3370 fsw. NMRI conducts hyperbaric 
research focused on the biomedical and physiological problems 
encountered by the fleet diver. 

3) Duke University's Hyperbaric Chamber - located in 
Durham, NC. This facility is capable of employing any desired 
gas at high pressures ranging from 150,000 ft of altitude to 
1000 fsw. Research 1s conducted on acoustic speech 
modification, high pressure nervous syndrome and carbon 
dioxide effects. 

A typical experimental hyperbaric facility is shown in figure 
ibe Joe 
Several factors must be considered when designing a 


hyperbaric chamber. NAVFAC DM-39 (1982) addresses the 





following design requirements: 

1) Pressure vessel design including material selection, 
wall thickness, welding requirements, connections, view ports, 
and doors. 

2) Fire protection and extinguishing systems. 

3) Electrical, utility, and communication systems. 

4) Life support systems including mixed gas breathing 
systems, carbon dioxide control, and oxygen replenishment. 

5) Temperature and humidity control. 

6) Power requirements. 

The objective of this program is to determine life support, 
temperature and humidity control, and power requirements. 
During most hyperbaric chamber designs, these factors are 
analyzed separately and compiled into a finished design. This 
program allows the user to input varying design considerations 
and provides a detailed output of design requirements. The 
most unique capability of this program is the use of a 
computer to determine the thermodynamic and psychrometric 
properties. Traditionally, psychrometric charts have been 
used to determine dehumidification, moisture removal, and 
ultimately energy requirements. Although specially developed 
psychrometric charts are available for the higher pressures 
required in a hyperbaric environment, they are time consuming 
and lack the accuracy of a computer program. Additionally, 
reading values off these charts requires non-linear 
interpolation. Therefore, the necessary relationships were 


developed to replace the psychrometric chart and simulate a 





hyperbaric chamber environment. 


INNER OUTER 
DECOMPRESSION DECOMPRESSION 
CHAMBER CHAMBER 


Divers Swimming and Working in Wetpot are 
Supported by Tenders in Igioo. 


Igloo is Primary Working, Living and Resting 
Space. ‘’Dry’’ Dives are also Conducted in 

this Space. 

Decompression Chamber is used for Sleeping 
and for Treatment of Decompression Sickness. 
Each Space in the Facility Can Be Pressurized 
Individually or in Combination With Other 
Spaces. 





Figure 1.1 Typical Experimental Hyperbaric Facility 
(Shilling et al. 1976) 


Objectives 


Based upon relevant aspects of human physiology, this 
program determines the ideal diving-gas mixture along with an 
acceptable range of diving-gas mixtures for a given depth or 
pressure. Additionally, the required ventilation rate without 
carbon dioxide absorbers or the amount of absorption material 
(sodasorb) for saturation diving is calculated. Various 
properties of the selected diving-gas mixture are determined, 
and by using thermodynamic properties of the gas, the heating 


and cooling requirements are calculated. Finally, the energy 


5 





requirements for dehumidification and reheating to a desired 
temperature are determined along with moisture removal 


requirements. 





CHAPTER II 


THEORY FOR DESIGN OF HYPERBARIC ENVIRONMENTS 


Basic Physiology 


A basic knowledge of human physiology is required to 
determine an acceptable diving-gas mixture for varying depths 
or pressures. While a given mixture may be acceptable at one 
pressure, it can cause discomfort or even be fatal at other 
pressures. The percentage of oxygen required to sustain life 
decreases aS pressure increases in order to maintain an 
acceptable oxygen partial pressure. 

The respiratory system of the body is responsible for 
oxygen uptake and carbon dioxide removal. Oxygen uptake is 
governed by two main factors listed by U. S. Navy Diving Gas 
Manual (1971): 

1. Respiratory Minute Volume 

2. Effectiveness in which the lung extracts oxygen from 

alveolar gas. 
Respiratory minute volume, which is defined as the product of 
the respiratory rate (breaths/min) and the tidal volume 
(volume/breath), determines the duration of exposure of blood 
to alveolar air. The effectiveness in which the lung extracts 
oxygen from alveolar gas depends upon 1) composition and 
volume of air supplied and 2) partial pressures. 

The transfer of gases takes place due to differences in 
partial pressures. The transfer of oxygen occurs with no 


complications if the oxygen partial pressure is between 0.16 
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and 1.6 atm for short periods of time. However, if the oxygen 
partial pressure is below 0.16 atm, anoxia (oxygen shortage) 
which may lead to fainting and possible brain damage can 
occur. Likewise, if the oxygen partial pressure exceeds 1.6 
atm, oxygen poisoning which may lead to convulsions and 
unconsciousness may occur. NOAA Diving Manual (1991) states 
that a 1.6 atm oxygen partial pressure can only be tolerated 
for up to 45 minutes and that the maximum oxygen partial 
pressure for long periods should not exceed 0.6 atm. NOAA 
Diving Manual (1991) also recommends that for long term 
exposures in a hyperbaric chamber or underwater habitat the 
oxygen partial pressure of the diving-gas should be between 
0.3 and 0.4 atm. Figure 2.1 shows the percentage of oxygen in 
the diving-gas mixture as a function of depth and oxygen 
partial pressure. Additionally, Figure 2.1 shows acceptable 
mixture ranges for a given depth and acceptable depth ranges 
for a given mixture. 

Under atmospheric conditions, the partial pressure of 
oxygen is 0.21 atm since the percentage of oxygen in air is 
approximately 21%. Therefore, at the surface, divers may 
safely breathe 100% oxygen for short periods (hours) since the 
oxygen partial pressure is only 1.0 atm, well below the 1.6 
atm upper limit. However, as depth increases, the pressure 
also increases and 1.6 atm can be quickly exceeded using air 
as the diving-gas. Equations 2.1 and 2.2 show the 
relationship between pressure (psia) and feet of sea water 


(fsw). 
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P(psia) = 


d(fsw) = 2.25P - 33.066 


where, 
33.066 fsw = 1 atm 
1 atm = 14.696 psia 
Additionally, the partial pressure of oxygen for given 


conditions can be determined from equations 2.3 and 2.4. 


. , a(fsw) 
EO 1 0) 80,) a 
_ , P(psia) 
PR Oe =) (wie S71) & 
eg 14.696 WER eyasih 


From equation 2.3 it can be shown that the greatest depth in 
which air can be safely used as a diving-gas without exceeding 
1.6 atm partial pressure is 218 fsw. 

In addition to quickly exceeding the allowable amount of 
oxygen, the presence of nitrogen renders air an undesirable 
breathing gas at depth. Under atmospheric conditions the 
partial pressure of nitrogen is low and no effects are felt. 
However, when the partial pressure of nitrogen is sufficiently 
high, nitrogen narcosis or rapture of the deep occurs. 


Nitrogen narcosis presents symptoms similar to alcohol 


se: 





intoxication and can be extremely dangerous while diving at 
great depths. Table 2-1 indicates the events which may occur 
while diving with high pressure nitrogen at depth (Randall 
1992). Therefore, air is not a suitable diving-gas and is not 
generally used for depths below 150 fsw. The advantage of 
using air when possible is that it is inexpensive, and a 
chamber can be ventilated continuously without the use of 


carbon dioxide absorbers which are addressed later. 


Table 2-1 Effects of Nitrogen Narcosis (Randall 1992) 


- 150 fsw light head, some euphoria 


- 200 fsw jovial, some dizziness 

- 250 fsw laughter may approach hysteria 
Simple tasks cannot be performed 
without mistakes 


depression 


may approach unconsciousness 





Since helium is an inert gas which eliminates the 
narcosis effect, it is often used at depths greater than 150 
to 200 fsw. However, helium does have the side effects which 
are listed in the NOAA Diving Manual (1991). One of these 
side effects is speech distortion which increases with depth. 
Speech distortion can be overcome by using helium speech 
scramblers. A second problem with helium is that its use 


results in rapid body heat loss since the thermal conductivity 
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of helium 1S six times that of air. To counter this problen, 
the diving-gas must be heated for dives below 800 fsw. A 
third problem associated with helium is the high pressure 
nervous syndrome (HPNS) which causes tremors in the hands and 
jerky movements of the limbs. This problem can be reduced by 
using a slow rate of compression. Alternatives to helium are 
also listed in the NOAA Diving Manual (1991). These 
alternatives include neon and hydrogen. Neon has’ the 
advantages over helium of a lower thermal conductivity and 
less speech distortion. However, neon is much more expensive 
and creates more breathing resistance. Hydrogen is easier to 
breathe at depth, but it causes more speech distortion than 
helium, has a higher thermal conductivity, and can be 
explosive. Other alternatives may be under development, but 
are beyond the scope of this report. 

This program recommends diving-gas mixtures based upon a 
helium-oxygen mixture. Smaller percentages of oxygen can be 
mixed with helium to provide a diving-gas mixture which can 
sustain life without the risk of oxygen poisoning at deeper 


depths. 


Gas Mixtures 
The ideal diving-gas mixture is determined from an oxygen 
partial pressure of 0.35 atm. Equation 2.5 calculates the 


percent oxygen in the ideal diving-gas mixture. 
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2 ideal P : 
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Since the diving-gases only contain helium and oxygen, 
equation 2.6 calculates the percent helium in the ideal 


diving-gas mixture. 


+He = 100 - %0, 


ideal 


The minimum and maximum allowable limits for oxygen percentage 
are calculated by equations 2.7 and 2.8 based upon oxygen 
partial pressures of 0.16 and 0.6 atm respectively as 


previously discussed. 


0.16 


atm 


%O = ( 


emin 


) x 100% 





- (0.6 


atm 


) x 100% 





If the maximum allowable oxygen percentage is greater 
ewan 21%, air is recommended as the diving-gas mixture. 
However, in most cases a hyperbaric chamber operates in an 
environment which necessitates the use of a helium-oxygen 
mixture. 


For a given diving-gas mixture, the pressure (pSia) and 


108) 





depth (fsw) ranges are calculated using equations 2.9 - 2.12. 








2.9 
_ ,16.0 
Pin = +0, ) (14.696) 
2.10 
60.0 
Po = (———) (14.696 
2 <M 
an (22:9) (33, 066) - 33.066 
0, 
22 





pew 


2 


) (33.066) - 33.066 


: 


These ranges are also shown in the previous Figure 2.1. 


Chamber Ventilation 

Ventilation is required to remove carbon dioxide and 
replenish oxygen to sustain life. When air is used, the 
required ventilation rate must be sufficient to flush the 
chamber to an acceptable level of carbon dioxide. This rate 
shall be referred to as the "free ventilation rate". The free 
ventilation rate is more than sufficient to provide enough 
oxygen to sustain life. Since air is relatively inexpensive, 
the free ventilation rate is desirable when utilizing air and 
carbon dioxide scrubbers are not required to keep the 


environment carbon dioxide free. 


14 





The free ventilation rate must be sufficient to maintain 
an environment at or below acceptable partial pressure levels 
of carbon dioxide. As with oxygen, the partial pressure of 
carbon dioxide increases with depth. For extended periods of 
time, the maximum allowable partial pressure of carbon dioxide 
should not exceed 0.005 atm. Figure 2.2 shows the maximum 
allowable CO, partial pressure with respect to mission 
duration. It is evident that as the mission duration 
increases, the smaller the carbon dioxide partial pressure 
must be in order to be tolerated. In the event that CO, 
partial pressure exceeds the allowable, hypercapnia or carbon 
dioxide excess can occur. Table 2-2 lists the four zones of 


hypercapnia. 


Table 2-2 Effects of Hypercapnia 


zone I No effect at low CO, partial pressures 


Zone II - Minor perceptive changes 
Zone III - Distracting discomfort 


Zone IV - Dizziness, stupor, unconsciousness 





Carbon dioxide production is determined by equation 2.13, 


CO, = (Rg) (Op, (Ng) 


ds. 
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where, 

C02, = carbon dioxide production (SCFM) 
Sag respiratory quotient (0.9) 

Oo. = oxygen consumption (SCFM) 


Cc 
Ns = number of divers 
Oxygen consumption is determined by the working conditions 


listed in Table 2-3. 


Table 2-3 Oxygen Consumption of Working Conditions 


Condition Description QO, consumption 
Resting 0.5 SLM / 0.0177 
Light Work 0.7 SLM / 0.0247 


Moderate Work 1.5 SLM 0.0530 


/ 
Heavy Work Den Omo Neen O07 OFF 
/ 


Severe Work 3.0 SLM 0.1060 





Figure 2.3 shows the oxygen consumption in Standard Liters per 
Minute (SLM) for each of the working conditions. Standard 
conditions exist at 14.696 psia and 70° F. 

In order to determine the free ventilation rate, a steady 
state condition and complete mixing are assumed. The 
in ft?/min is 


volumetric flow rate of air at depth (V,;,-) 


determined by equation 2.14. 
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where, 
m = the mass flow rate of CO, (lbgo7/minute) 


R 


the gas constant for CO, (35.11 lb/lbmole-R) 

T = the initial temperature (R) 

PCO, = the partial pressure of CO, (psf) 

(0.005 atm or 10.58 psf is recommended for times 
greater than 8 hours) 


The ventilation rate (S.:_) in Standard Cubic Feet per Minute 


air 


(SCFM) is calculated by equation 2.15. 


GS. = ( GVarr) Geers) ) 


aa 14.696 


The mass flow rate of CO, is calculated by equation 2.16. 


: ; “le 
(14.696 psia) (144 psia/psf) (CO, ) i 
(35.11 ft-lb/1lbmole-R) (530 R) 
where, 14.696 psia and 530 R are standard conditions. 
Equation 2.16 can be simplified to equation 2.17. 
Zena, 


m = (0.1137) (CO, ) 


ibs) 





where m is expressed in lbgo2/minute. 

Since a steady-state condition is assumed, it is 
important to determine how long it takes to reach this 
condition at the beginning of a mission or when conditions 
change. The time required to reach a 95% steady-state 


condition can be derived from equation 2.18. 


- (Ve) (t) ) 


PCO, = (ee) [1=e we] 2.18 


alr 





where V,, = the interior volume of the chamber ee) 


Rearranging to solve for t yields equation 2.19. 


Ma 


alr 


PCO, , 
MRT 2.19 


se (A) In [a= | 


air 


By substituting in the values for PCO, (0.005 atm), 
R (35.11 ft-lb/lbmole-R), and a 95% steady-state condition, 


equation 2.19 becomes equation 2.20. 


mos) oes) (14.7) (144) (V,,.) . 


eS Veh = 
eee ‘F ia’ (35.11) (mT) GAG 


air 


Carrying out the multiplication yields equation 2.21. 
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=O Oley jet Vv ; 
a. ) (Vem) ) (an py — ¢ £0+2863) (Vase) 
Vair mT Za 


where the factor 0.0167 converts t from minutes to hours. 


Carbon Dioxide Removal 

In most cases, the working depth or pressure requires a 
helium-oxygen mixture for reasons previously stated. In these 
cases scrubbers are used to remove carbon dioxide from the 
Chamber. Therefore, pure oxygen is injected but care must be 
taken to ensure proper mixing with helium in the chamber. 
This mixing occurs while the gas mixture is recirculating 
through the dehumidification and scrubbing system. 

The amount of pure oxygen required to sustain the divers 
(@5-) in Standard Cubic Feet per Minute (SCFM) is determined 


from equation 2.22. 


0, = (0, ) (Ny) 


where oxygen consumption (0,.) was previously determined from 
the diver working conditions. The amount of CO, produced 


(lbgo9) in pounds per minute is calculated by equation 2.23. 


1beo, = (m) (ty) (60) 
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where t, 1s the mission duration and the factor of 60 converts 
minutes to hours. 

NAVFAC DM-39 (1982) lists three chemicals used as carbon 
dioxide absorbents. These chemicals and their properties are 
shown on Table 2-4. Baralyme is no longer used in hyperbaric 
chambers and Lithium Hydroxide is too expensive and difficult 
to handle for routine use. Therefore, Sodasorb is the 
chemical absorbent used in this program. The amount of 
Sodasorb (LB,) required for each mission is calculated from 


equation 2.24 assuming a 50% efficiency. 


LD¢o, 
ee 0.245 


Although it is possible to obtain higher absorbent 
efficiencies, canister design, temperature, humidity, and 
moisture condensation must be considered. The 50% efficiency 
is a conservative estimate used for design purposes. During 
actual chamber operations, chemical detection tubes must be 
used to accurately determine when to change carbon dioxide 
absorbent canisters. The number of carbon dioxide absorbent 
Canisters required for a mission will depend on _ the 


Capabilities of individual hyperbaric chambers. 


Diving-Gas Mixture Properties 


Properties of the diving-gas mixture must be accurately 


determined in order to calculate energy and dehumidification 
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Table 2-4 Characteristics of Carbon Dioxide Absorbents 
(NAVFAC DM-39 1982) 


Absorbent 


Characteristic Baralyme Lithium Sodasorb 
ee Sy Groxide — 


Density (lb/ft?) 65.4 28.0 


Theoretical Co, 0.39 0.92 
absorption (lb CO,/1b) 


Useful CO, absorption 0.195 -46 
(based on 50% efficiency) 
(lb CO,/1b) 





requirements which will be discussed later. Properties of 
concern include molecular weight, specific heat, and density 
of the selected diving-gas mixture. 

The molecular weight of a substance is defined as the 
weight in pounds of one pound-mole of the substance. The 
molecular weight of the diving-gas mixtures is dependent upon 
the percentages of oxygen and helium in the mixture and 
therefore, decreases as the percent oxygen in the mixture 
decreases. This in turn decreases the density which minimizes 
the energy requirements. If the diving-gas is air, the 
molecular weight is 29.0 lb/lb-mole. Otherwise, the mole 
fractions of oxygen and helium in the diving-gas mixture must 
first be determined. Mole fractions (X) are determined by 
equations 2.25 and 2.26. The sum of the mole fractions must 


always equal 1. 
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29 Xo, 


Once mole fractions are determined, the molecular weight in 


lb/lb-mole can be calculated from equation 2.27. 


MW = [(X_) (MWo)] + [(Xy.) (My) | 


where, 
MWo2 = 32.00 1lb/lb-mole 
MWie = 4.0026 1lb/lb-mole 


The specific heat at constant pressure (C in units got 


oue 


BTU/lb-mole, of the diving-gas mixture must be known for 


determining energy requirements. Cp is defined by equation 
Brea . 
2.28 
@ =C 
Pp OT Pp 


where, h is the enthalpy and T is temperature. 
Given the mole fractions previously determined, Cy in BTU/1lb- 


mole can be calculated from equation 2.29. 


eee) 91 (%,,.) (C,, )] 


24 





6.95 BTU/1lb-mole 


Q 
ll 


4.965 BTU/lb-mole 
When air is used as the diving-gas, Cy is 6.95 BTU/l1b-mole. 
Table 2-5 lists the molecular weights and specific heats of 


common helium/oxygen mixtures and air. 


Table 2-5 Properties of Common Diving-Gas Mixtures 


Diving-Gas Mixture © MW 


$05/%He) (BTU/1—mole) (1b/lb-mole) 
1/99 4.985 4.283 


5/95 5.064 5.403 
10/90 .164 6.802 
15/85 ws 8.202 


20/80 .362 9.602 


40/60 .759 15.202 


Air -950 -000 





Density of air is calculated using the ideal gas equation 
of state since air will only be used down to 150 fsw (81.36 


psia). The ideal equation of state is shown as equation 2.30. 


PV = nR,T 
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where, 
P = absolute pressure (psf) 


V 


volume (cu ft) 


n number of moles 
R, = Universal Gas Constant (1544.0 ft-1lb/1bmole-R) 
T = absolute temperature (R) 


Equation 2.30 can be simplified and rearranged to solve for 


density in equation 2.31. 


_ (MW) (P) 


Density of helium-oxygen mixtures is determined from the 
virial equation of state. Dymond and Smith (1980) present the 


virial equation of state shown as equation 2.32. 











where, the values of P,R,, and T are the same as those used 
above, and 
Vm = molar volume (cu ft/1lb-mole) 


1 = the first virial coefficient 


B(T) = the second virial coefficient 


law the third virial coefficient 


D(T) the fourth virial coefficient 
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It is common practice to calculate the density of gases based 
upon only the first three virial coefficients since data 
beyond this range are scarce. Dymond and Smith (1980) state 
that uncertainties in the second virial coefficient are 
relatively small (usually less than 1%) but are much larger 
for the third virial coefficient. Additionally, the third 
virial coefficient is much more temperature dependent. 
Dymond and Smith (1980) developed the most recent table 
for determining the second and third virial coefficients for 
gas mixtures. Equation 2.33 calculates the second virial 


coefficient of a gas mixture. 


Peewee 5. Ween) + B, (1) (X,) (xX) + B,, (T) (x,)* 


where, 
B,(T) = the second virial coefficient of the mixture 
X, = the mole fraction of the first pure gas 


the mole fraction of the second pure gas 


x 
NS 
ll 


= the second virial coefficient of the first pure gas 


109) 
amy 
ra 

| 


B55 = the second virial coefficient of the second pure gas 
B15 = the mixed virial coefficient 
Equation 2.34 calculates the third virial coefficient of a gas 


mixture. 


2d 





eee Ca AG) eG. 1) (Xx, ° (X,) 
fmt eres () 4X5.) (X,)'* +:C55 (7) 4X.) 2 2.34 


where, 
C,(T) = the third virial coefficient of the mixture 
Cii1 = the third virial coefficient of the first pure gas 
Cii2 = the contribution of two molecules of the first pure 
gas and one molecule of the second pure gas 
Cioo = the contribution of one molecule of the first pure 
gas and two molecules of the second pure gas 
Due to the absence of data for determining the third virial 
coefficients for helium-oxygen mixtures, the density is 
calculated based upon the first two virial coefficients only. 
Tables 2-6 lists the virial coefficients for oxygen, heliun, 
and a helium-oxygen mixture respectively from Dymond and Smith 
(1980). Since this is the most recent data available, it is 
obvious that very little data exists to determine the mixed 
virial coefficient for helium-oxygen mixtures. 
By substituting the values from Tables 2-6 into equation 
2.33, the second virial coefficient in SI units (cm?/g-mole) 


is calculated from equation 2.35. 


B= ((11.44) (X,.)7] - [(8.8) (Xy.) (Xo) 1 
Se lG 24 ee! 2235 
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Table 2-6 Virial Coefficients of Helium, Oxygen, and 
Heliox Mixtures 
(Dymond and Smith 1980) 


hacia leGoertlcrents 


Helium Oxygen Mixture 
B S B ‘8 


10.76 ~12764 
Od ‘ 1230 
11.44 ‘ 1163 


als Bese «7; ° 1146 


B= cm?/g-mole = cm® /g-mole* 





In order to use this coefficient, it must be converted into 
the English units of ft?/lb-mole. This conversion is made by 


equation 2.36. 


at oe (454 J) 
28316.85cm? Ib 2.36 


(0.0160) (B,-) 


Solving the virial equation of state for molar volume (ft?/1b- 


mole) yields equation 2.37. 
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where K is a constant defined as R,T/P. 


Finally, density is calculated by equation 2.38. 


SE 


where MW was calculated in equation 2.27 and molar volume was 
calculated in equation 2.37. 

The U. S. Navy Diving Gas Manual (1971) lists the 
densities of standard diving-gas mixtures at varying 


conditions. The resulting error from these equations is less 


than 0.2%. 


Psychrometric Charts 

Psychrometric charts are used to compute heating and 
cooling requirements along with dehumidification and moisture 
removal requirements. Standard psychrometric charts can be 
used for calculating the effects of heating and cooling moist 
air at atmospheric conditions but cannot be used with diving- 
gas mixtures at much greater depths. Although these charts 
are available for pressures up to 600 psia (approx. 1300 fsw) 
in the U. S. Navy Diving Gas Manual (1971), they are time 
ensuing and lack the accuracy of this computer program. 
Therefore, it is necessary to develop relationships which can 
be used in computer simulation. Additionally, reading values 


off these charts requires non-linear interpolation between 100 


psia curves. Figure 2.4 shows the psychrometric chart in the 
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U. S. Navy Diving Gas Manual (1971) developed for gas mixtures 
at pressures from 100 to 600 psia. 

Given the dry bulb temperature and saturation lines at 
constant pressure, the moisture mole fraction (moles water 
vapor/moles dry gas) can be determined. The dry bulb 
temperature is defined by ASHRAE (1967) as the Fahrenheit 
temperature of moist gas at rest with respect to the 
temperature-measuring element. The saturation lines at 
constant pressure represent 100% relative humidity at that 
pressure. 

The moisture mole fraction in moles water vapor/mole dry 
gas is determined by equation 2.39 if the pressure of the 


environment and saturation pressure are known. 


where, 
X35 = the mole fraction for water vapor at saturation 
(moles water vapor/mole dry gas) 


P, = the vapor pressure of water in saturated moist air 


(psia) 
P = the absolute pressure of the chamber (psia) 
The values of P, have been tabulated by ASHRAE (1967) for 
various temperatures. These values are reproduced in Table 


2-7. Equation 2.40 is a best fit curve through these points 
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Table 2-7 Vapor Pressure of Water in Saturated Moist Air (psia) 
ASHRAE (1967) 


Absolute Pressure (Ps) 


Temp (F) Ps (psia) Temp (F) Ps (psia) 
4Q OF iio 72 OeSs6G06 
41 0.12646 ive 0.40190 
42 Omir oes 74 0.41564 
43 Oph GeO es 0.42979 
44 0.14194 as OU4aaa 52 
49 0.14746 TM ORAS S 
46 Oivosled 78 | DAs 
4°7 Chole) SNC ire 79 0.49066 
48 Ol Goi, 80 Os5 070! 
49 0.17148 81 OR eyaete 2 
ale) Ol77 2 eve 0.904112 
owl, Oh ey4ere 3) So OVDecge2 
me Oh Shilens, 84 Oe iar 2 
a0 0.19883 85 0.59604 
o4 0.20630 86 0.61540 
= 5) lees 1 See Ee O)6)635)-0,.0 
o6 SF aye ilers: 88 ORG 8 fo 
a OFZ 3006 89 0.67678 
o8 0.23849 90 0.69838 
og Or 70 91 OF7 2059 
60 Oh. alee: Je 0) aeak Ss abG 
61 O 26040 25) OL elsicial 
62 Ope 2 94 0.790901 
Se! 0.28488 99d 0.81564 
64 O20 500 96 OTo4 Ue 
oS, 30504 a4 OPSCralel 
66 0.31636 98 0.89388 
67 Oro (a0 215) QO: Ozer 
68 0.33900 100 0.94959 
69 0.35064 
70 0.86304 
al ORS TAD LOM 
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determined from PCMATLAB (1989) with a correlation coefficient 


oi, O.9999965. 


Peel 285 x 1017) T° + (9.8683 x 1071") 7 


mcs slos x 10-7) T° + (1.9916 x 10°) 2.40 
+ (1.0488 x 10°7)T + (2.2080 x 1072) 


where T = the chamber temperature (° F). 

The U. S. Navy Diving Gas Manual (1971) states that for 
total mixture pressures up to 100 psia, no correction for the 
effect of total pressure on the water-vapor saturation 
pressure is_ used. However, above 100 psSia, Poynting's 
equation iS used as an estimated correction for total 
pressure. Glasstone (1947) derives Poynting's Equation shown 


aS equation 2.41. 


_ Vv, (AP) (P,) 


AP. 
RT 2.41 


where V, is the specific volume of the saturated liquid. 
Values of V, are tabulated by ASHRAE (1967) and are reproduced 
in Table 2-8 in units of (ft°?/lb). To be used in Poynting's 
Equation, V, must be expressed in (ft?/mole). This is easily 
accomplished by multiplying V, by the molecular weight of 
water (18.0158 lb/lb-mole). Then, the corrected saturation 


pressure will be calculated from equation 2.42. 
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Table 2-8 specific Volume of Saturated Liquid 


Temp (F) 


40 
Al 
Ae 
43 
bt 
495 
46 
4°7 
48 
49 
oO 
Ol 
OK 
Od 
o4 
ore) 
O6 
i 
08 
o°9 
60 
Cul: 
oe 
63 
64 
69 
66 
67 
68 
69 
70 
TN 


Vy 


7 


0.01602 
0.01602 
Ore re0< 
CeO NGUZ 
0.01602 
0.01602 
HO leO2 
0.016Ce 
0.01602 
0.01602 
Ge0T6Ge2 
GrG1602 
0.01602 
0.01603 
0.01603 
© e1603 
Gre1503 
0.01603 
0.01603 
0.01603 
0.01603 
0.01604 
0.01604 
0.01604 
0.01604 
0.01604 
0.01604 
0.01619 
0.01609 
0.01609 
0.01605 
0.01609 


ASHRAE (1967) 


Sweet roe, Ol Uiml 
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Temp (F) 


Te 
1 
74 
79 
76 
Cl 
78 
Ws) 
80 
81 
ee 
83 
84 
85 
86 
87 
88 
89 
ZO 
G1 
ve 
os 
94 
99 
wo 
97 
98 
eke 
100 


yi 
0.01606 
0.01606 
0.01606 
0.01606 
0.01606 
CEOS Oi 
OO 1607 
GEO moe 
OO ke Oy 
0.01608 
0.01608 
0.01608 
O-O1606 
0.01609 
0.01609 
SONOS, 
OCS 1G 
0.01610 
C20 lies 
Oh Ge ING: 
OO Tor 
OPO Tes IL IE 
DOV r | 
GO OHIS Lee 
Sroore 
O20 LG 12 
0.01612 
0.01613 
ORO Ne os 





where P,,; 1s the uncorrected saturation pressure. The 
moisture mole fraction of water vapor at saturation (X,,) is 


then calculated from equation 2.39. 
The molar enthalpy of moisture content at saturation (H) 


in (BTU/mole dry gas) is calculated from equation 2.43. 


H = (20,000) X,, 


Enthalpy is a thermodynamic property defined by equation 2.44. 


igh (10M caw <a 


where, 
U = internal energy 
P = absolute pressure 


V = volume 
Determination of enthalpy is important for calculating energy 


requirements. 


Dehumidification and Energy Requirements 


In order to determine dehumidification requirements, the 


initial and desired relative humidities must be known. Given 
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the initial and desired relative humidities, the moisture mole 
fraction and molar enthalpy of moisture content can be 


determined for both conditions by equations 2.45 - 2.48. 


X,, = (RH;) (X,,) 2.45 
X,, = (RH) (X,5) 2.46 
H, = (RH;) (#) 2.47 

2.48 


H, = (RH) (H) 


where RH; and RH, are the initial and desired relative 
humidities expressed in fractions, and X,, and H were 
previously determined. Once the moisture mole fraction at the 
desired relative humidity is known, the maximum coil 
temperature required for dehumidification can be calculated. 


Equation 2.49 calculates the new saturation pressure (P,). 


Cx) <P) 
eel Xe 2.49 
2 


CW, 





Then, the maximum coil temperature (T,) in (°F) can be 
determined as the intersection of saturation pressure curve at 


feo Or from equation 2.50. 


ieee. 4135 x 10°) P - (8.4538 x 10°) BP,’ 
ele 6 35x10") P.> miS. 1910 x10*) P.2 250 
mo 407 x 107) P+ (7.6106) 


Equation 2.50 was calculated using PCMATLAB (1989) and has a 
correlation coefficient of 0.999974. 

Once the coil temperature is determined, the moist gas 
mixture enthalpies at the initial and coil temperatures are 
Calculated. Equations 2.51 and 2.52 calculate the moist gas 
mixture enthalpies in (BTU/1b) at the initial and coil 


temperatures respectively. 


2.51 
= LG) (T;) + i) 
% MW 

2.52 
wg L(G) (Te) +H] 
2 MW 


where the specific heat and molecular weight were previously 
determined as properties of the diving-gas mixture. 

The energy (Q4) in (BTU/min) required to dehumidify the 
chamber to the desired relative humidity is calculated from 


Saadci1on 2.53. 
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where, V.;, 1s either the free ventilation rate for air 
systems or the recirculation flowrate for helium-oxygen 
systems. Since the coil temperature is generally colder than 
the desired | Gdiar Meee mixture, the diving-gas must be 
reheated to a desired temperature. The change in enthalpy 
Puameche Coll to desired temperature (Tg) is calculated by 


equation 2.54. 


meee rer, nee 


Then, the energy required to reheat the diving-gas (Q,) to the 


desired temperature is calculated by equation 2.55. 


QO, — p VairAh 


The total energy (Q,) required by the system _ for 


dehumidification and reheat is calculated by equation 2.56. 


Cas Oat 0, 


However, energy lost to heat transfer through the chamber 
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walls depends on insulation and is not considered in this 


analysis. It would be an additional energy requirement. 


Moisture Removal Rate 

In all dehumidification processes, moisture removal must 
be addressed. Due to the rapid cooling of the moist air, 
water vapor condenses and must be removed from the chamber. 
In order to determine the moisture removal rate, the humidity 
ratios at the initial and desired relative humidities must be 
determined. Equations 2.57 and 2.58 calculate the humidity 
ratios in pounds water vapor per pounds dry air (1b,/lbgy) at 


the initial (W,) and desired (W,) relative humidities 


respectively. 
sos 
en) (18.0158) ; 
W, SS ee 
MW 

Pa iB ye) 

(ee) 37,0158) 

W, — 


MW 


where 18.0158 is the molecular weight of water. The moisture 
removal rate (MRR) in (lb/day) can then calculated by equation 


2.09. 


a 


MRR = (W, - W,)V,;,P (1440)  lb/day 2.59 


where, V,;, is either the free ventilation rate or the 
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recirculation flowrate and the factor of 1440 converts 
lb/minute to lb/day. This moisture removal rate results from 
the dehumidification process only. The moisture produced by 
the carbon dioxide absorbents is not included in this 


analysis. 
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CHAPTER III 


DESCRIPTION OF SIMULATION 


Simulation Structure 

This program is designed to simulate hyperbaric chamber 
environments such as decompression chambers, submarine lockout 
and pilot epembeney diving bells, underwater habitats, and 
ocean simulation facilities under different conditions. The 
theory and equations developed in Chapter II are incorporated 
into this interactive computer program which is written in 
Fortran 77 for a 386 personal computer. In order to design a 
hyperbaric chamber with this program, the input listed in 
Table 3-1 1s required. Given this input, the program 
calculates values and provides the output listed in Table 3-2. 

A flowchart for this program is included as Figure 3.1 


and the source code is provided as Appendix B. 


Subroutine Description 


This program makes use of several subroutines’ to 
calculate the required quantities. A description of each 


subroutine follows in the paragraphs below. 


Subroutine "CONVERT" converts depth to pressure. For all 
subsequent calculations, this program uses absolute pressure. 
If the user chooses to input a working depth in feet of 
Seawater (fsw), "CONVERT" converts depth to pressure using 


equation 2.1. This value is returned to the main program and 
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Papte 3-1 Required Program Input 


Requirement Units 


1) Working Depth fsw (feet of seawater) 
or Pressure psia 


2) Initial Temperature oF 
3) Desired Temperature uk 
4) Initial Relative Humidity % 

5) Desired Relative Humidity % 

6) Number of Divers N/A 
7) Internal Volume of Chamber ees 
8) Diver Working Conditions N/A 
9)* Mission Duration hrs 


10)“ Recirculation Flowrate ft?/min 


* - These values are only required if CO, scrubbers are 
going to be used during the mission. 





used for all future calculations. 


Subroutine "GASMIX" calculates the minimum and maximum oxygen 
percentages for the pressure of interest based upon oxygen 
partial pressures of 0.16 and 0.6 atm respectively from 
equations 2.7 and 2.8. If the maximum acceptable oxygen 
percentage is greater than 21%, the ideal diving-gas mixture 
is air. All further calculations are then based upon air as 
the diving-gas mixture. If the maximum acceptable oxygen 


percentage is less than 21%, the program determines an 


acceptable oxygen percentage using an IF-THEN statement. 
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Mable 3=2 Program Output 


Output Units 
1) Diving-Gas Mixture sHe/%05 


2) Oxygen Required SCFM 
or Ventilation Rate 


3") Sodasorb Required lbs 
4) Time to Reach 95% Steady-State hrs 


5) Molecular Weight lb/1lb-mole 


6) Gas Mixture Density Nb/tt- 


7) Maximum Coil Temperature ae 

8) Dehumidification Energy Requirements BTU/min 
9) Reheat Energy Requirements BTU/min 
10) Moisture Removal Rate lb/day 


- This value is only provided if CO, scrubbers are 
used during the mission. 





Subroutine "'GASRANGE" calculates a pressure and depth range 
based upon equations 2.9 - 2.12 after an acceptable diving-gas 
mixture has been selected. This information is useful if the 
same diving-gas mixture is used for multiple depths. If the 
user is planning to operate outside these limits, a different 


diving-gas mixture should be selected. 


Subroutine "FREERATE" determines the ventilation rate required 
to keep the carbon dioxide partial pressure at an acceptable 


level when air is used as the diving-gas. First, "FREERATE" 
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calculates the carbon dioxide production based upon the 
working conditions and equation 2.13. In order to maintain a 
CO, partial pressure below 0.005 atm, the volumetric flow 
rates of air at depth and at the surface are calculated from 


equations 2.14 and 2.15. 


Subroutine "ABSRATE" calculates the amount of oxygen required 
to sustain life, the amount of CO, produced, and the amount of 
Sodasorb required to maintain an acceptable environment when 
the diving-gas is a helium-oxygen mixture. The amount of 
oxygen required is calculated by equation 2.22 while the 
amount of CO, produced in pounds is calculated by equation 
2.23 and the amount of Sodasorb required is calculated by 
equation 2.24. The user may wish to design different CO, 


scrubbing systems based upon the amount of CO, produced. 


Subroutine "STSTATE" calculates the time required to reach a 
95% steady-state condition. Since steady-state conditions are 
assumed, it is important to determine the time required to 
reach this condition. The time required to reach the 95% 
steady-state condition is calculated by equation 2.21. This 
information is useful when beginning a mission or changing 


conditions. 


Subroutines "MOLWT", "SPECHEAT", and "DENSITY" calculate the 
properties of the selected diving-gas mixture. If the diving- 


gas mixture is air, the molecular weight and specific heat are 
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constant values of 29.0 lb/lb-mole and 6.95 BTU/lb-mole 
respectively. 

However, for helium-oxygen mixtures, "MOLWT" calculates 
the molecular weight from equation 2.27 and "SPECHEAT" 
calculates the specific heat from equation 2.29. This 
information is used in subsequent calculations. 

"DENSITY" calculates the density of the diving-gas for 
both air and helium-oxygen mixtures. Since air is only used 


at low pressures, "DENSITY" determines the density of air from 


the ideal equation of state, equation 2.31. When helium- 
oxygen mixtures are used, "DENSITY" calculates density from 
the virial equation of state, equation 2.32. The second 


virial coefficient for the mixture is determined from equation 
2.35 and converted to English units in equation 2.36. MThis 
value is then used to determine the molar volume from equation 
2.37, and finally the density is calculated by equation 2.38. 


This value is used in subsequent energy calculations. 


Subroutine "PSYCHART" calculates the moisture mole fraction 
(X,,) and the enthalpy of moisture content (H) for the diving- 
gas mixture. First, the vapor pressure of water in saturated 
moist air (P,) is calculated from equation 2.40. If the 
working pressure is less than 100 psia, no correction is 
required. However, for pressures above 100 psia, P, is 
corrected using Poynting's equation 2.41 where the specific 
volume of the saturated liquid (V,) is read off of Table 2-8 


using an IF-THEN statement. The corrected value of P, is then 
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calculated by equation 2.42. The values of X,, and H are 


calculated from equations 2.39 and 2.43 respectively. 


Subroutine "XWSH" calculates the initial and desired moisture 
mole fractions (X,, and X,.) and molar enthalpies (H, and H>). 
These properties are determined by the initial and desired 
relative humidities. The values are calculated based upon X,,, 


and H by equations 2.45 - 2.48. 


Subroutine "COILTEMP" calculates the maximum coil temperature 
required to dehumidify the chamber to the desired relative 
humidity. The moisture mole fraction at the desired relative 
humidity (X,5.) was calculated by equation 2.46. The new 
saturation pressure (P,) is determined by solving equation 
2.39 for P, which leads to equation 2.49. Given P,, the coil 


temperature (T,) is calculated by equation 2.50. 


Subroutine "MGMENTH" calculates the moist gas mixture 
enthalpies entering the cooling coil (H,) and leaving the 
cooling coil (H,) by equations 2.51 and 2.52 respectively. In 
this program, the temperature entering the cooling coil is 
assumed to be the initial temperature (T;) and the temperature 
leaving the cooling coil is assumed to be the maximum coil 
temperature (T,).- These values may differ in an actual 


cooling system. 


Subroutine "ENERGY" calculates the energy required for 
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dehumidification and reheat to the desired temperature. The 
energy required for dehumidification (Q,) is calculated by 
equation 2.53. If the diving-gas is air, the free ventilation 
rate (V,;,) calculated by "FREERATE" is used. For helium- 
oxygen mixtures, the recirculation flowrate input by the user 
is used as V,;,- In order to determine the energy required to 
reheat to a desired temperature (Q,), the change in enthalpy 
from the coil temperature to a desired temperature must be 
calculated by equation 2.54. Then, Q, is calculated in a 
Similar manner as Q, by equation 2.55. The total energy 


required (Q,) is calculated by equation 2.56. 


Subroutine "MOISTREM" calculates the moisture removal rate 
required for the above dehumidification. First, the humidity 
ratios at the initial (W,) and desired (W,) relative 
humidities are calculated by equations 2.57 and 2.58. Then, 
the moisture removal rate is calculated by equation 2.59 where 
again, V.;, is the free ventilation rate for air or the 
recirculation flowrate for helium-oxygen mixtures. The 
calculations performed by these subroutines are incorporated 
into the main program, and the results are displayed to the 


user as program "output". Three examples are discussed in 


detail in the next chapter. 
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CHAPTER IV 


APPLICATION OF SIMULATION 


This program can be used not only to simulate a 
hyperbaric chamber but also a submarine diver lock-out chamber 
or an underwater habitat. The following examples demonstrate 
this program's usefulness for determining the requirements of 


these cases. 


Example 1: Hyperbaric Chamber 


This example simulates a hyperbaric chamber used for a 
deep dive. As a model, the U. S. Navy's Ocean Simulation 
Facility at the Naval Coastal Systems Center in Panama City, 
Fl was be used. This facility as shown in Figure 4.1 is 
capable of pressurization down to 2250 fsw and is described in 
detail by Shilling et al. (1976) and U. S. Navy Diving Manual 
(1991). The input conditions for this example are given in 


Table 4-1. The calculation steps are demonstrated below: 


1) Since a working depth was input, it was converted to 


a pressure. By equation 2.1, 


(2000 + 33.066) ) 


14.696) 
33.066 


P= ( 


v 
rT 


903.6 pSia 
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Figure 4.1 Ocean Simulation Facility 
(Shilling et al. 1976) 


2) The ideal gas mixture for this depth was caiculated 


by equations 2.5 and 2.6. 


tO, = (03.6 )” (100%) 


14.696 


$0, = 0.57% 


%He 
%He 


HOO. = 0.57 
99.43% 
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Table 4-1 Hyperbaric Chamber Simulation Input 


Depth (fsw): 2000 
Initial/Desired Temperatures (° F): 90/80 
Initial/Desired Relative Humidities (%): 90/50 


Number of Divers: 6 


Internal Chamber Volume (ft3): 3300 


Recirculation Flowrate (cfm): 100 
Mission Duration (hrs): 240 


Diver Working Conditions: Light Work 





Therefore, the ideal gas mixture was determined to be 


0.57% O, and 99.43% He. 


3h) The minimum and maximum oxygen percentages were 


calculated by equations 2.7 and 2.8. 


Oe 6 


=e ee Coe 
tO, ae ) 
14.696 
$0, = 0.26% 
0.6 
Se a ee Fekete 
SO, ‘7903.6, ) 
14.696 
$0, = 0.98% 
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Since only mixtures within these levels are acceptable, 
air could not be used and a 0.5% 0, and 99.5% He diving- 
gas mixture was selected. Therefore, carbon dioxide 
scrubbers were used and a computation of free ventilation 


Was not needed. 


4) The pressure and depth ranges for the gas mixture 





were calculated from equations 2.9 - 2.12. 
+ Gaal 6 
Paes ok a= ) (14.696) 
Pun = 470.3 psia 
ee Se 
Pee = (14.696) 
Jno, ee dag) sh 5 Saegs) lel 





aaa = Le * =) (335066)8-7 33.066] 


1025.0 fsw 


Q. 


mim 





a= petals: (33.066) - 33.066] 


O52 
nna 


3934.9 fsw 


5) Since a helium-oxygen mixture was used with carbon 
dioxide scrubbers, the amount of pure oxygen injected 


into the chamber was calculated by equation 2.22. 
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Oo. = 0.7 SLM or 0.0247 SCFM for light work 


O, = (0.0247) (6) 
O, = 0.1482 SCFM 
6) The amount of carbon dioxide produced and the 


scrubber requirements were calculated by equations 2.17 


emee2ze2s — 2.24. 


(0.1137) (059) (0.1482) 
0.01517 Ibe. /min 


= 
T 


= 
T 


Ibeo, = (0.0157) (240) (60) 


1beo, = 218.4 Ib 


27 .218.4 
1Ds ‘9. 245) 
De 389 1, 4 
7) The time required to reach 95% steady-state was 


determined from equation 2.21 where V,;, was determined 


fiom 2.14. 
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» © « GE aap @&@ —= =p 


a! 
-_ a 


EO vOrS17) (3551 lye49 or 





ae 
a 10.58 
Vee ee) eben 
£ = (-0.0167) (3299) inpy-—(0:2863) (27.7) _ | 
D7 7 (0.0152) (549.67) 
t = 5.96 Ars 


8) Next, the diving-gas properties were determined. 


Molecular weight was calculated from equations 2.25-2.27. 


x = 0.5 — 0.005 
2 00 
9935 7 
Xe = “100° = 0.995 


(emOOSee2 OO) (0; 99s) 0026) 4 
4.143 lb/lb-mole 


s 3 


Specific heat was calculated from equation 2.29. 
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= [(0.005) (6.95) + (0.995) (4.965) ] 
= 4.975 BTU/lb-mole 


Density was calculated from equations 2.35 - 2.38. 


[(11.44) (0.995)? - (8.8) (0.995) (0.005) 
- (16.24) (0.005) 2] 


Oy 
l 


11.2187 cm?/g-mole 


by 
il 


CE OnrG 0) atnlie 2, Galleys) 


Oy 
i 


0.1805 ft?/lb-mole 


Oy 
T 


_ (1544) (549.67) 
(903.6) (144) 


K=6.5226 


60 2260 6 226) (0.1 805 
2 
0.6185 ft?/lb-mole 


ES 
ll 


p = 4-143 
6.6984 


0.6185 lb/ft? 


ca 
ll 
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9) The psychrometric properties were calculated from 
equations 2.39 - 2.48. The psychrometric chart created 


for this case is presented as Figure 4.2. 


P, = 0.69838 psia (from Table 2-7) 
eee eo-efolo cu ft/lb (from Table 2=8) 
V, = 0.2880 cu ft/mole 


Ap. = [(0.69838) (0.2880) (903.6) (144)] 
: [(1544) (549.67) ] 


AP 0.0308 psia 


0.69838 + 0.0308 
0.7292 psia 


On Zou. 
a SOS 605 Oe 207 


x 
i 


0.00081 moles water vapor/mole dry gas 


WS 


(X,, is shown as point A on Figure 4.2) 


(0.9) (0.00081) 


T 


wl 


= 0.00072 moles water vapor/mole dry gas 


< 


wl 
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Moisture Mole Fraction 


1000 0 


aS 


<000 0 
£000 0 
VOOO 0 
S000 0 
9000 0 
£000 0 


8000 0 
6000 0 


OTO0 0 
TROO 0 
e100 0 


&TO0 0 


(moles water vapor/mole dry g 


VLOO'O 


S1LOO 0 


(MSJ) OOO? 4OJ JAeYD OlMeuIOAYyoOASG Sp eins 


(1) oanyestodutay, qung Aad 


: | UoIJBoTIprUNnYyaq : | , | : 
eee 72 -BuNe09 | EO OED ll (ot ae ee 





08 





(X,, 1S shown as point B on Figure 4.2) 


x 
| 


PeOrs) (0.00087) 


< 
U 


w2 = 9.00041 moles water vapor/mole dry gas 


(X,5 1S shown as point D on Figure 4.2) 


H = (20,000) (0.00081) 
H = 16.2 BTU/mole dry gas 


ES (ORM 
H, = 14.6 BTU/mole dry gas 
He —{0.5) (1652) 
H, = 8.1 BTU/mole dry gas 


10) The maximum coil temperature was determined from the 
psychrometric chart as the intersection of X,. and the 
100% saturation curve. This point is at approximately 
70° F and is shown as point C on Figure 4.2. 


= O 
ih Oa 13 
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lls) The energy requirements were calculated from 


equations 2.51 - 2.56 


Ryo 309 0) ahaa 
‘ 4.143 


Hldie,(S Jou Jie 


ce 
Tl 


2 (4 Se 7 0) + 8.0) 
: 4.143 


86.0 BTU/1b 


ce 
T 


Og (Onoheay4LOO) All 6"— 86.0) 


O, = 1583.4 BTU/min 


a ESET GO SO) 
eis 


Ah = 12.0 BTU/1b 
Cer= Onoda) LOG) (lz .0) 
Op -942.2 BLU min 


O, = 1583.4 + 742.2 


Or 


2325.6 BTU/min 
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12) The moisture removal rate was calculated from 


equations 2.57 - 2.59. 


18.0159 
a COMOOI0g 2))..( ————__— — 
- ( ) ¢ 4.143 
Wee O,00sise1D /1b, 

18.0158 
(ae OM OO O41) (= = 
2 ( ) 4.143 
W, = 0.00178 1b,/1b, 


MRR = (0.00313 - 0.00178) (100) (0.6185) (1440) 
20), 2 b/ day 


. 


The actual program printout is provided as Table 4-2. The 
slight discrepancies between the example and the actual output 
were caused by round-off error obtained while working the 


example. 
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Table 4-2 Output for Hyperbaric Chamber Simulation 


HYPERBARIC CHAMBER SIMULATION 


Input Conditions: 
Pressure (psia): 


Initial/Desired Temperature (° F): 


Initial/Desired Relative Humidity (%): 


# of Divers/Working Conditions: 
Chamber Volume (ft?): 
Flow Rate (ft?/min): 


Mission Duration (hrs): 


Output Conditions: 

Diving-Gas Mixture (%He/%0,) : 
Oxygen Required (SCFM): 

CO, Produced (lbs): 

Sodasorb Required (lbs): 

Time to 95% Steady-State (hrs): 
Molecular Weight (1b/l1b-mole): 

Gas Mixture Density (1b/ft?): 
Maximum Coil Temperature (° F): 
Dehumidification Energy (BTU/min) : 
Reheat Energy (BTU/min): 

Total Energy Requirement (BTU/min): 


Moisture Removal Rate (lb/day): 
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903.6 


90.0/80.0 


90.0/50.0 
6.0/2 
3300.0 
100.0 


240.0 


99.5/0.5 
0.148 
218.7 
892.7 
6.0 
4.143 
0.6184 
70.2 
1563.5 
724.3 
2287 ec 


$256.2 





Example 2: Submarine Diver Lock-out Chamber 


This example simulates a diver lock-out chamber aboard a 
submarine used to decompress divers after a mission. The 
input conditions for this example are given in Table 4-3. 


Table 4-3 Diver Lock-out Chamber Simulation Input 


Depth (fsw): 
Initial/Desired Temperatures (° F) 
Initial/Desired Relative Humidities (%): 


Number of Divers: 


Internal Chamber Volume (ft): 


Recirculation Flowrate (cfm): 10 
Mission Duration (hrs): 4 


Diver Working Conditions: Resting 





The calculation steps are demonstrated below: 
1) Since a working depth was input, it was converted to 


a pressure. By equation 2.1, 


= ((400 + 33.066) ) | 


14.696) 
33.066 


U 
T 


192.5 psia 


2) The ideal gas mixture for this depth was calculated 


by equations 2.5 and 2.6. 
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%0, = (53-5) (100%) 


14.696 


%O, = 2.67% 


ll 
Nw 


%*He = 100 - 2.67 
%He = 97.33% 


Therefore, the ideal gas mixture was determined to be 
2.67% O5 and 97.334 He; 


3) The minimum and maximum oxygen percentages were 


calculated by equations 2.7 and 2.8. 


0, | _ ( 0.16 
ae JS PAYS: ) 
14.696 


1.2% 


1 (ALIORS 


2min 


7 0.6 
$0, = ee, 


14.696 


$0, = 4.58% 


Since only mixtures within these levels air acceptable, 
air could not be used and a 3% O, and 97% He mixture was 


selected. Therefore, carbon dioxide scrubbers were used 
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and a computation of free ventilation rate was not 
needed. 


4) The pressure and depth ranges for the gas mixture 








were calculated from equations 2.9 - 2.12. 
P= | ae (14.696) 
Pin = 78.4 psia 
60.0 
ae = A 50 ) (14.696) 
Jog AE Se foKepte! 





ote Re eeroes) ro SOG: | 


320 
aan 


143.3 fsw 





(cA an a a (Se q066 y- 33.066] 


ax 


628.2 fsw 


5) Since a helium-oxygen mixture was used with carbon 
dioxide scrubbers, the amount of pure oxygen injected 
into the chamber was calculated by equation 2.22. 


Oo, = 0.5 SLM or 0.0177 SCFM 
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G5 (0. 0177) (2) 
O, = 0.0354 SCFM 
6) The amount of carbon dioxide produced and the 


Scrubber requirements were calculated by equations 2.17 


aWemee.2S5—- 2.24. 


CO. 137 ).( 0. 99200 -.03 SA) 


mM — 
m= 0.0036 Lbeo /min 
1beo, = (0.0036) (4) (60) 
Ibe, = 0.86 Ib 
0.86 
1b. = (2a. 
z ‘9.245? 
iene ae aime le, 
7) The time required to reach 95% steady-state was 


determined from equation 2.21 where V.,. was determined 


from 2.14. 
vy. = (0.0036) (35.11) (544.67) ] 
air 
Oso 
V Soa chee fact amata) 
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WOmz2ec65 ) (6551) ] 


100 c 
ea) ail (0.0036) (544.67) 





bE =P On O67) 


t = 0.77 Ars 


8) Next, the diving-gas properties were determined. 


Molecular weight was calculated from equations 2.25 


Ze! « 


2 ~=400 





97.0 _ 
Sn 2 ay DOSY 


MW = [(0.03) (32.00) + (0.97) (4.0026) J 
MW = 4.843 lb/lb-mole 


Specific heat was calculated from equation 2.29. 


eo 03) (6.95) 4 (0.97) (4.965) ] 


C) 
T 


5.025 BTU/lb-mole 


CY 
T 


Density was calculated from equations 2.35 - 2.38. 
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elie) 500. 977 ) = — 8, 8)(0.97)40'..03 ) 
=e howe 4) (0.03) 7%] 


by 
T 


B = 10.4932 cm?/g-mole 


B= (0.0160) (10.492) 
B= 0.1679 ft?/lb-mole 
_ (leap aa) 

(192.5) (144) 
K = 30.34 


_ 30.34 + 730.542 + 4(30.54) (0.1679) 


V,, = 
2 
Var 3045s. / LD-mole 





Je 4.843 
30.51 
oe O)e Ab svi eeyoy/imia 


9) The psychrometric properties were calculated from 
equations 2.39 - 2.48. The psychrometric chart created 
for this case is presented as Figure 4.3. 

P. = 0.59604 psia (from Table 2-7) 

Pea omoroOo secu ft/lb (from Table 2-8) 


Wiaelsceceo Cu £t/mole 
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Moisture Mole Fraction 


Ol O07 O6 O8 OL O9 OG 
O00 0 
| OO © Se oe eae a se ees anne n a'wneg ce -eeeee see ata as Meee ns oan wa seams anne Ne etna sede ae eum eaetene ce ae2- ov eaaas donee anesa sens a-matemeeineay 22 Se swans and motion ten eeneesin ae aeeeavenia 25 asada va-wansuusue sp saec caer tahoe d= sei eee ee ee aes can des 
yeoyoey 9 
TAG OS. road oe: eee e renee eae I oc Sopa epee + Seo so reeecrere Re ecco ce Se eee d ne, oie sopctora eopere eee eee YY ee ee Reece era haere mens Peer [eeweene s+ MEMES ARnea tense erawesseene 


S010 78 


(moles water vapor/mole dry gas) 


"OOD 


ENNSiig eek Aeele) tev n2)16 pg ware) ob 8 Kom Ab rey Csr ele amc a alre Ny 
({) eunyetodtsay, qing Aaq 


oT] eorIprunyeaq 
; 9» SUI[OOD 


OV 
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Ap. = [(0-59604) (0.2896) (192.5) (144)] 
; [ (1544) (544.67) ] 


AP 0.00569 psia 


= 0.59604 + 0.00569 


0.60173 psia 


Oe oUl > 
ne EZ oP 100017 3 


< 
T 


0.00314 moles water vapor/mole dry gas 


WS 


(X,, 1S shown as point A on Figure 4.3) 


be 
I 


1 (Oro (0700314) 


Ww 


be 
Il 


“4 = 0.00282 moles water vapor/mole dry gas 


(X,,; 1S shown as point B on Figure 4.3) 


ora ONG EO OOS 4) 


x 
T 


v2 = 0.00188 moles water vapor/mole dry gas 


70 





(X,5 1S shown as point D on Figure 4.3) 


H = (20,000) (0.00314) 
= 62.71 BTU/mole dry gas 


a 


iso = 0) 0) NG aay al), 


H, = 56.44 BTU/mole dry gas 
ah, (lias) (52 al) 
H, = 37.63 BTU/mole dry gas 


10) The maximum coil temperature was determined from the 
psychrometric chart as the intersection of X,. and the 
100% saturation curve. This point is at approximately 
70° F and is shown as point C on Figure 4.3. 

T, = 70° F 


11) The energy requirements were calculated from 


equations 2.51 2656: 


Det>. 025)7(085) 7% 56.44] 
2 4.843 


99.85 BTU/1b 


ty 
Tl 


aL 





S025 (7 0)eeteed. 6 3 | 
4.843 


H. = 80.40 BTU/I1b 


Ie ms 


OF — ie 15 Sr eos oes se> 80.40) 
QO, = 30.9 BTU/min 
i eelicrOZ ease 70) 
4.843 
Ah = 8.3 BTU/I1b 
Oe = HOe eset tno) C8 3) 
QO, = 13.2 BTU/min 
OP Oe St oer 
QO, = 44.1 BTU/min 
2) The moisture removal rate was calculated from 


equations 2.57 - 2.59. 


18.0158 
Wa = OOO eo 
a ( ) { 4.843 ) 
W, = 0.01049 1b,/1b, 


a2 





138.0156 
W, = (0.00188) ( ———— 
ne 188) ( Tae 
W, = 0.00699 1b,/1b, 


MRR = (0.01049 — 0.00699) (10) (0.1587) (1440) 
MRR = 8.0 lb/day 
The actual program printout is provided as Table 4-4. The 


slight discrepancies between the example and the actual output 
were caused by round-off error obtained while working the 


example. 
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Table 4-4 Output for Diver Lock-out Chamber Simulation 


DIVER LOCK-OUT CHAMBER SIMULATION 


Input Conditions: 


Pressure (psia): 192.5 
Initial/Desired Temperature (° F): 85.0/78.0 
Initial/Desired Relative Humidity (%): 90.0/60.0 
# of Divers/Working Conditions: 250) fal 

Chamber Volume (ft): 100.0 
Flow Rate (ft?/min): TOs 0 
Mission Duration (hrs): 4.0 


Output Conditions: 


Diving-Gas Moisture (%He/%0,) : S707 3°30 
Oxygen Required (SCFM): 0.035 
CO, Produced (lbs): 0.9 
Sodasorb Required (lbs): 355 
Time to 95% Steady-State (hrs): 0.8 
Molecular Weight (1b/lb-mole): 4.843 
Gas Mixture Density (lb/ft?): 0.1587 
Maximum Coil Temperature (° F): 70.0 
Dehumidification Energy (BTU/min): 30.9 
Reheat Energy (BTU/min): Secor 
Total Energy Requirement (BTU/min): 44.1 
Moisture Removal Rate (lb/day): 8.0 
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Example 3: Underwater Habitat 


This example simulates a shallow underwater habitat. The 
French Conshelf II is used as a model. Conshelf II, shown as 
Figure 4.4, was a shallow water habitat operating in the Red 
Sea in September of 1963. The breathing gas was air and the 
mission duration was 1 month. Miller and Koblick (1984) 
G@escribe this mission in detail. The input conditions for 
this example are given in Table 4-5. The calculation steps 


are demonstrated below: 


Table 4-5 Underwater Habitat Simulation Input 


Depth (fsw): 50 
Initial/Desired Temperatures (° F): 95/80 
Initial/Desired Relative Humidities (%): 95/50 


Number of Divers: 4 


Internal Chamber Volume (ft?): 2000 


Diver Working Conditions: Moderate Work 





1) Since a working depth was input, it was converted to 


a pressure. By equation 2.1, 


= (90 + 33.066) , | 


14.696) 
33.066 


P = 36.9 psia 
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Figure 4.4 Conshelf II Habitat, Starfish House 
(Miller and Koblick 1984) 


2) The ideal gas mixture for this depth is calculated by 


equations 2.5 and 2.6. 


Os: 
J629 
14.696 


13.9% 


(100%) 


| 
“~—, 


$0, = 


$0, 


HOOe= 13:59 
86.1% 


%He 
%He 
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3) The minimum and maximum oxygen percentages were 


calculated by equations 2.7 and 2.8. 


0.16 
%0 =o) (ceo (100% 
2min 36.9 s 
14.696 
$0, = 6.37% 
0.60 
20. = (00S 
2max a 36.9 a ) 
14.696 
$0, = 23.9% 


Since any mixture within these levels is acceptable, air 
fle O>)ewas selected as the diving-gas. Therefore, a 
free ventilation rate was used to keep the chamber free 


of carbon dioxide and CO2 absorbers were not required. 


4) The free ventilation rate at depth and from the 


surface was calculated from equations 2.14 - 2.17. 
Oo, = 1.5 SLM or 0.0530 SCFM for moderate work 
Core — 00 2) ehOr 0.550) (4) 


0.1908 lb/min 


CO, 
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= 
I 


(Om 7) (061908) 
0.0217 1bgo /min 


3 
T 


Pom Oz vena Same) (529.67 ) | 


Vi 
a 10.58 
Viir = 39-9 ft?/min 
og = 1699.9) (36.9) (529.67) | 
ae [(14.696) (554.67) ] 
Se = oo eS eel 
5) The time required to reach a 95% steady-state 


conditions was determined from equation 2.21. 


a 2000 noe, coos) (39.9) 
E ( 0.0167) (35-5) inf1 7 SE EGY. 
tc =2.5 Ars 


6) Next, the diving-gas properties were determined. 

Since the diving-gas was air, the molecular weights and 
specific heats were 29.0 lb/lb-mole and 6.95 BTU/1b-mole 
respectively. Density was calculated from the ideal 


equation of state, equation 2.31. 


vo 





_ [(29.0) (36.9) (144) ] 


7 ea) SS Aaa ))) 
p = 0.18002 Ib/ft? 
7) The psychrometric properties were calculated from 


equations 2.39 - 2.48. The psychrometric chart created 
for this case is presented as Figure 4.5. 

P, = 0.81564 psia (from Table 2-7) 
Since the operating pressure was less than 100 psia, no 


correction was required for the saturation vapor 


pressure. 


0.81564 
OF BO. = OB ESo4 


0.0226 moles water vapor/mole dry gas 


< 
Tl 


WS 


(X,, 1S shown as point A on Figure 4.5) 


SOS 5) 400226 ) 


wl 


0.0215 moles water vapor/mole dry gas 


x 
T 


wl 


(X,, 1s shown as point B on Figure 4.5) 


TES, 





Moisture Mole Fraction 
(moles water vapor/mole dry gas) 


(MSJ) OG 1IOJ JIeYD OliyOWIOIUOASG Gy ainsiy 


(J) FInyetodue] qing Aid 


OU OO | Q6 O08 oe OF, OS Ov 
OOO 


qeoyey 


L0°0 
_ uolqeotjiptiunyag 
@ TSUT[OOD 
Z0'O 
c0'0 
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care (On50) (Gn02 ogy 


0.0113 moles water vapor/mole dry gas 


x 
T 


w2 


(X,5 1S shown as point D on Figure 4.5) 


H = (20,000) (0.0226) 
H = 452.0 BTU/mole dry gas 


H, = (0.95) (452.0) 
H, = 429.4 BTU/mole dry gas 
H, = (0.50) (452.0) 
H, = 226.0 BTU/mole dry gas 


8) The maximum coil temperature was determined from the 
psychrometric chart as the intersection of X,, and the 
100% saturation curve. This point is at approximately 
73.7° F and is shown as point C on Figure 4.5. 

| ae <a al 


9) The energy requirements were calculated from 


equations 2.51 - 2.56. 
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Wher 2S) (See One) 
: 29.0 


Bi au aly 


a 
T 


= 43,95) (Veo) 8 PAs) 
= 29.0 


25.5 870/ ib 


ee 
Tl 


QO, = (0.18002) (39.9) (37.6 - 25.5) 
QO, = 86.9 BTU/min 


ee eee rT) 
29.0 


Ah = 1.51 BTU/1b 


QO, = (0.18002) (39.9) (1.51) 
QO. = 10.8 BTU/min 


QO. = 86.9 + 10.8 
0, 


97.7 BTU/min 


10) The moisture removal rate was calculated from 


equations 2.57 - 2.59. 
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SLSTE LOSS 
W, = (0.18002) (—————— 
a=" ) ( Sar0 
W, = 0.01336 1b,/1b, 

Leh, alse! 
W, = (0.18002) (—————— 
aa ) ( SG ) 
W, = 0.00702 1b,/1b, 


Vere Oeols3s6 = @.00702) (39.9) (0.18002) (1440) 
MRR = 65.5 lb/day 
The actual program printout is provided as Table 4-6. The 


slight discrepancies between the example and the actual output 
were caused by round-off error obtained while working the 


example. 
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Table 4-6 Output for Underwater Habitat 


UNDERWATER HABITAT SIMULATION 


Input Conditions: 


Pressure (psSia): 36.9 
Initial/Desired Temperature (° F): 95.0/80.0 
Initial/Desired Relative Humidity (%): 95.0/50.0 
# of Divers/Working Conditions: 4.0/3 
Chamber Volume (ft?): 2000.0 


Output Conditions: 


Diving-Gas Mixture: AIR 
Free Ventilation Rate at Depth (cfm): 39.9 
Free Ventilation Rate (SCFM): 95.8 
Time to 95% Steady-State (hrs): 7S, 
Molecular Weight (1b/l1b-mole): 29.0 
Gas Mixture Density (lb/ft?): 0.18002 
Maximum Coil Temperature (° F): 73.8 
Dehumidification Energy (BTU/min): 87.0 
Reheat Energy (BTU/min): 10.8 
Total Energy Requirement (BTU/min): 97.8 
Moisture Removal Rate (lb/day): 65.4 
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CHAPTER V 


SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS 


A FORTRAN computer program has been developed to simulate 
a hyperbaric chamber environment. Although hyperbaric 
chambers have been used for many years, increasing complexity 
requires innovative design procedures. This program allows 
the user to design life support systems, temperature and 
humidity control equipment, and determine energy requirements 
for hyperbaric chambers, submarine diver lockout chambers, and 
underwater habitats. By using this program, the design 
engineer can Simulate a variety of scenarios and determine the 


design criteria for a specific chamber. 


Conclusions: 

1) Acceptable diving-gas mixtures are obtained by ensuring 
the oxygen partial pressure of the mixture remains between 
0.16 and 1.6 atm for short periods of time and below 0.6 atm 
for extended operations. 

2) In order to eliminate the narcosis effects of nitrogen 
on the diver, helium-oxygen mixtures are used below 150 fsw. 

3) The partial pressure of carbon dioxide (CO,) must be 
maintained below 0.005 atm during extended missions in order 
to eliminate discomfort and possible injury to the divers. 
This can be accomplished by ventilation of the chamber or by 
using a CO, absorption system. Ventilation is'~ only 


recommended when air is used as the diving-gas due to the cost 
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of helium. 

4) A CO, absorption system consists of canisters filled 
with sodasorb. The diving-gas mixture flows through these 
canisters and the CO, is removed. Although this program can 
be used to estimate the duration of the CO, canisters if the 
amount of sodasorb contained is known, a chemical detection 
system should always be used to determine when to replace the 
absorbent material. 

5) The molecular weight and specific heat at constant 
pressure of the diving-gas mixture can be determined from the 
component mole fractions. 

6) The density of helium-oxygen mixtures is determined from 
the virial equation of state. Three virial coefficients are 
required to accurately calculate density. Although, the first 
two virial coefficients have been accurately determined for 
helium and oxygen, only one second virial coefficient for 
helium-oxygen mixtures has been determined. Additionally, no 
third virial coefficients for helium-oxygen mixtures have been 
determined to date. Therefore, extremely accurate density 
calculations for helium-oxygen mixtures are not obtainable at 
this time. 

7) Psychrometric charts have been used for years to 
determine dehumidification and cooling requirements. 
Relationships were developed to compute moisture mole 
fractions (X,,) allowing a computer simulation without using 
the psychrometric charts. This program is easier to use and 


more accurate than psychrometric charts. 
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8) Enthalpy of moisture content (H) is easily calculated 
once X,, is determined. Enthalpy is an important 
thermodynamic property for calculating energy requirements. 

9) Dehumidification and energy requirements can be 
determined for varying conditions if the initial and desired 
temperatures and relative humidities are known. However, the 
total energy requirement neglects heat loss due to heat 
transfer between the chamber and_~ the environment. 
Additionally, the coil temperature required to perform the 
desired dehumidification and the moisture removal rate can be 


determined. 


Recommendations: 

1) Pressure vessel design of a hyperbaric chamber could be 
incorporated into a computer program. NAVFAC DM-39 (1982) 
provides guidance for pressure vessel design of hyperbaric 
chambers and the Pressure Vessel Handbook (1977) simplifies 
the calculations of pressure vessel design for desired 
conditions. 

2) Fire protection and extinguishing systems along with 
other safety design criteria could be incorporated into a 
computer program. Additionally, electrical, utility, and 
communication systems could be considered. NAVFAC DM-39 
(1982) provides guidance for these systems. 

3) Different diving-gas mixtures could be investigated with 
additional facility requirements incorporated into’ the 


program. The NOAA Diving Manual (1991) lists some 
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alternatives to helium and research into the effects and 
advantages of these gases is ongoing. 

4) In order to accurately calculate the density of diving- 
gas mixtures, the second and third virial coefficients need to 
be accurately determined. Although these coefficients exist 
for helium and oxygen, they have not been accurately 
determined for the mixtures. Texas A & M University has the 
facilities to accurately determine these coefficients if 
funding is available. 

5) Heat transfer between the hyperbaric chamber and the 
environment is not considered in this program. A detailed 
study of materials and insulation would be necessary prior to 


incorporating it into a computer program. 
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APPENDIX A 


NOMENCLATURE 


B - second virial coefficient, ft?/mole 

C - third virial coefficient, ft®/mole? 

C02, - carbon dioxide production, SCFM 

Cy 7 specific heat at constant pressure, BTU/1lb-mole 

D - fourth virial coefficient, ft? /mole? 

a - working depth, feet of seawater (fsw) 

dnax ~ maximum allowable depth for a diving-gas mixture, fsw 

dain ~ Minimum allowable depth for a diving-gas mixture, fsw 

H - molar enthalpy of moisture content at saturation, 

BTU/mole dry gas 

- moist-gas-mixture enthalpy at initial temperature, BTU/1b 

- moist-gas-mixture enthalpy at coil temperature, BTU/1b 

H, - molar enthalpy of moisture content at initial conditions, 
BTU/mole dry gas 

H, - molar enthalpy of moisture content at desired conditions, 
BTU/mole dry gas 

lb.g9 - Carbon dioxide produced, lb 

lb, - sodasorb required for absorption, lb 

m - mass flowrate of carbon dioxide, 1b/min 

MRR - moisture removal rate, lb/day 

MW - molecular weight of diving-gas mixture, lb/lb-mole 

n - number of moles 

N,; - number of divers 

2 


¢ 7 oxygen consumption, SCFM 
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oe emygen required to sustain divers, SCFM 


P - absolute pressure, psia 


Prax ~ Maximum allowable pressure for a diving-gas mixture, 
psia 

Pain ~ minimum allowable pressure for a diving-gas mixture, 
psia 


PP5> - partial pressure of oxygen, atm 


PCO, - partial pressure of carbon dioxide, atm 


fe] — Vapor pressure of water in saturated moist air, psia 
HE; geal ~ Percentage of helium in ideal diving-gas mixture 


S05;Ggeal ~ Percentage of oxygen in ideal diving-gas mixture 

SOonax ~ Maximum acceptable oxygen percentage 

S0onin ~ Minimum acceptable oxygen percentage 

Qyg - energy required for dehumidification, BTU/min 

Sees energy required to reheat from coil to desired 
temperature, BTU/min 

eee EOtal energy required, BITU/min 

R - gas constant for CO,, 35.11 lb/lbmole-R 

RH, - desired relative humidity, % 

RH; - initial relative humidity, % 


Ra 7 respiratory quotient 
R, ~ universal gas constant, 1544 ft-lb/lbmole-R 


S - ventilation rate of air at standard conditions, SCFM 


air 
T - initial temperature, ° F 
t - time to reach 95% steady-state conditions, hr 


T. - maximum coil temperature for dehumidification, ° F 


Ty - desired temperature, ° F 
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t, - mission duration, hr 
U - internal energy, BTU 
V - volume, ft? 


V - volumetric flowrate at depth, ft?/min 


air 


Ven ~ volume of the chamber, ft? 


Cc 


V, - specific volune, ft?/1b 


V.. ~ molar volume, ft?/mole 


= 
_ 
1 


humidity ratio at initial relative humidity, 1b,/lb, 
W, - humidity ratio at desired relative humidity, 1b,,/1b,, 
Xue ~ mole fraction of helium 
Xoo ~ mole fraction of oxygen 
Xg ~ mole fraction of water at saturation, 
moles water vapor/mole dry gas 
X17 7 moisture mole fraction at initial conditions, 
moles water vapor/mole dry gas 
X29 - moisture mole fraction at desired conditions, 


moles water vapor/mole dry gas 
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APPENDIX B 


SOURCE CODE 
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aa 


DQAAAANAINANAAAAAANANAANAAANNNAANQAANANDNAANRAAANDAANAARANRARAANDN 


MAJOR REPORT - "A COMPUTER SIMULATION OF A HYPERBARIC 
CHAMBER ENVIRONMENT." 


REAL VAR, DEPTH, PRESS, TFI, TRI, TFD, RHI, RHD, 
+NUMDIV, CHVOL, FLOWRATE, DUR, VAR1, O2CONS, PERO2, 
+PERHE, O2MIN, O2MAX, DMIN, DMAX, PSIAMIN, PSIAMAX, VAIR, 
+SURAIR, O2REQ, TIME, MW, CP, DENS, BC, BF, XWS, XWl, 
+XW2, H, H1, H2, HA, HC, CT, QD, QR, QT, MRR, AIR, LBCO2, 
+LBS, CASES 


VAR = A VARIABLE TO DETERMINE WHETHER THE USER WISHES TO 
CALCULATE THE GAS MIXTURE BASED UPON DEPTH OR 
PRESSURE. 

DEPTH = A WORKING DEPTH USED TO CALCULATE THE DIVING-GAS 

MIXTURE IN FEET. 
PRESS = PRESSURE AT DEPTH IN PSIA. 
TFI = THE INITIAL TEMPERATURE OF THE CHAMBER IN DEGREES 


FARENHEIT. 

TRI = THE INITIAL TEMPERATURE OF THE CHAMBER IN DEGREES 
RANKINE. | 

TFD = THE DESIRED TEMPERATURE OF THE CHAMBER IN DEGREES 
FARENHEIT. 

RHI = THE INITIAL RELATIVE HUMIDITY IN THE CHAMBER. 

RHD = THE DESIRED RELATIVE HUMIDITY IN THE CHAMBER. 


NUMDIV = NUMBER OF DIVERS IN THE CHAMBER. 
CHVOL = THE INTERNAL VOLUME OF THE CHAMBER IN CUBIC 
PERT 
FLOWRATE = THE FLOWRATE INSIDE THE CHAMBER USED FOR 
RECIRCULATING THE DIVING-GAS MIXTURE WHEN CO2 
SCRUBBERS ARE USED. 
DUR = THE MISSION DURATION IN HOURS. 
VAR1 = A VARIABLE USED TO DETERMINE THE WORKING 
CONDITIONS. 
O2CONS = OXYGEN CONSUMPTION IN CUBIC FEET PER MINUTE. 
PERO2 = THE PERCENT OF OXYGEN IN THE DIVING-GAS MIXTURE. 
PERHE THE PERCENT OF HELIUM IN THE DIVING-GAS MIXTURE. 
O2MIN = MINIMUM ACCEPTABLE OXYGEN PERCENTAGE IN THE 
DIVING-GAS. 
O2MAX = MAXIMUM ACCEPTABLE OXYGEN PERCENTAGE IN THE 
DIVING-GAS. 
DMIN = THE MINIMUM DEPTH IN FEET FOR WHICH A GIVEN 
DIVING-GAS CAN BE USED. 
DMAX = THE MAXIMUM DEPTH IN FEET FOR WHICH A GIVEN 
DIVING-GAS CAN BE USED. 
PSIAMIN = THE MINIMUM PRESSURE IN PSIA FOR WHICH A GIVEN 
DIVING-GAS CAN BE USED. 
PSIAMAX = THE MAXIMUM PRESSURE IN PSIA FOR WHICH A GIVEN 
DIVING-GAS CAN BE USED. 
VAIR = VENTILATION RATE AT DEPTH IN CUBIC FEET PER 
MINUTE. 
SURAIR = VENTILATION RATE AT THE SURFACE IN STANDARD 
CUBIC FEET PER MINUTE. 
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OOOO Ora aadgmyraaandanananantonrn Yo” (wermnO087OM 4200000 


QO 


O2REQ = OXYGEN REQUIRED TO SUSTAIN DIVERS IN STANDARD 
CUBIC FEET PER MINUTE IF CO2 ABSORBERS ARE USED. 
TIME = THE TIME REQUIRED FOR THE CHAMBER TO REACH 99% 
STEADY-STATE CONDITIONS. 
MW = THE MOLECULAR WEIGHT OF THE DIVING-GAS MIXTURE 
(LB/LB-MOLE) . 
CP = THE SPECIFIC HEAT OF THE DIVING-GAS MIXTURE 
(BTU/MOLE) . 
DENS = THE DENSITY OF THE DIVING-GAS MIXTURE (LB/CU FT). 
BC = THE SECOND VIRIAL COEFFICIENT (CU CM/G-MOLE) . 
BF = THE SECOND VIRIAL COEFFICIENT (CU FT/LB-MOLE) . 
XWS = THE MOISTURE MOLE FRACTION OF WATER VAPOR AT 
SATURATION (MOLES WATER VAPOR/MOLE DRY GAS). 
XW1 = THE MOISTURE MOLE FRACTION AT THE INITIAL RELATIVE 
HUMIDITY. 
XW2 = THE MOISTURE MOLE FRACTION AT THE DESIRED RELATIVE 
HUMIDITY. 
H = THE MOLAR ENTHALPY OF MOISTURE CONTENT AT SATURATION 
(BTU/MOLE DRY GAS). 
Hl = THE MOLAR ENTHALPY OF MOISTURE CONTENT AT THE 
INITIAL R.H. 
H2 = THE MOLAR ENTHALPY OF MOISTURE CONTENT AT THE 
DESIRED R.H. 
HA = THE MOIST-GAS-MIXTURE ENTHALPY AT INITIAL 
CONDITIONS (BTU/LB DRY GAS). 
HC = THE MOIST-GAS-MIXTURE ENTHALPY AT DESIRED 
CONDITIONS. 
CT = THE COIL TEMPERATURE REQUIRED FOR DEHUMIDIFICATION 
IN DEGREES FARENHEIT. 
QD = THE ENERGY REQUIRED FOR DEHUMIDIFICATION (BTU/MIN). 
QR = THE ENERGY REQUIRED FOR REHEAT (BTU/MIN). 
QT = THE TOTAL ENERGY REQUIRED FOR DEHUMIDIFICATION AND 
REHEAT (BTU/MIN). 
MRR = THE REQUIRED MOISTURE REMOVAL RATE FOR THE 
DEHUMIDIFICATION PROCESS (LB/DAY). 
LBCO2 = AMOUNT OF CO2 PRODUCED (LBS). 
AIR = A VARIABLE TO DETERMINE WHETHER OR NOT AIR IS 
BEING USED AS THE DIVING-GAS. 
LBS = AMOUNT OF SODASORB REQUIRED FOR CO2 ABSORPTION 
(LBS) . 
CASES = THE NUMBER OF CASES TO BE CALCULATED. 


THIS PORTION OF THE PROGRAM WILL ALLOW THE USER TO INPUT 
THE NECESSARY PARAMETERS FROM A DATA FILE. 


OPEN (UNIT=3,FILE='VAR.DAT!,STATUS='OLD') 
OPEN (UNIT=4, FILE='OUTPUT.DAT! , STATUS='UNKNOWN' ) 


READ (3,*) CASES 


DO 500 I = 1,CASES 
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READ (3,*) VAR, DEPTH, TFI, TFD, RHI, RHD 
READ (3,*) NUMDIV, CHVOL, FLOWRATE, DUR, VARI 


IF (VAR.EQ.0) THEN 
CALL CONVERT(DEPTH, PRESS) 
GO TO 10 

ELSE 
PRESS=DEPTH 
GO TO 10 

ENDIF 

TRI=TFI+459.67 


IF (VAR1.EQ.1) THEN 
O2CONS=0.5/28.3 
ELSEIF (VAR1.EQ.2) THEN 
O2CONS=0.7/28.3 
ELSEIF (VAR1.EQ.3) THEN 
O2CONS=1.5/28.3 
ELSEIF (VAR1.EQ.4) THEN 
O2CONS=2.0/28.3 
ELSE 
O2CONS=3.0/28.5 
ENDIF 


THIS PORTION OF THE PROGRAM CALCULATES THE IDEAL AND 
ACCEPTABLE DIVING-GAS MIXTURES. 


CALL GASMIX (PRESS,O2MIN,O2MAX) 
IF (O2MAX.GE.21.0) THEN 


ATR=1 
PERO2=21.0 


ELSE 
AIR=0 

IF (O2MAX.GE.15.0) THEN 
PERO2=15.0 

ELSEIF (O2MAX.GE.10.0) THEN 
PERO2=10.0 

ELSEIF (O2MAX.GE.5.0) THEN 
PERO2=5.0 


ELSEIF (O2MAX.GE.3.0) THEN 
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PERO2=3.0 
ELSEIF (O2MAX.GE.1.0) THEN 
PERO2=1.0 
ELSEIF (O2MAX.GE.0.5) THEN 
PERO2=0.5 
ELSE 
PERO2=0.1 
ENDIF 
PERHE=100.0-PERO2 
CALL GASRANGE(PERO2,DMIN, DMAX, PSIAMIN, PSIAMAX) 


ENDIF 


THIS PORTION OF THE PROGRAM WILL CALCULATE THE REQUIRED 
VENTILATION RATE WITH OR WITHOUT CARBON DIOXIDE 
ABSORBERS. 
IF (AIR.EQ.1) THEN 

CALL FREERATE(02CONS,NUMDIV, TRI, PRESS, VAIR, SURAIR) 
ELSE 

CALL ABSRATE(NUMDIV,O2REQ,02CONS, DUR, LBCO2, LBS) 
ENDIF 
CALL STSTATE(CHVOL, O2CONS,NUMDIV, TRI, TIME) 


THIS PORTION OF THIS PROGRAM CALCULATES THE PROPERTIES 
OF THE DIVING-GAS MIXTURE. 


CALL MOLWT(PERO2,MW, AIR) 
CALL SPECHEAT(PERO2,CP,AIR) 
CALL DENSITY (PERO2,TRI, PRESS, MW, DENS, BC, BF, AIR) 


THIS PORTION OF THE PROGRAM DETERMINES THE PSYCHROMETRIC 
PROPERTIES. 


CALL PSYCHART(TFI,TRI,PRESS, XWS,H) 
CALL XWSH(RHI,RHD,XWS, XW1, XW2,H,H1,H2) 
CALL COILTEMP (XW2, PRESS, CT) 


CALL MGMENTH (TFI,CT,CP,H1,H2,MW,HA,HC) 


97 





Qa 


a2 


15 


20 


25 


30 


35 


40 


45 


50 


DS 


60 


THIS PORTION OF THE PROGRAM CALCULATES THE ENERGY 
REQUIREMENTS FOR DEHUMIDIFICATION AND REHEAT ALONG WITH 
MOISTURE REMOVAL REQUIREMENTS. 


CALL ENERGY (DENS, AIR, VAIR, FLOWRATE,HA,HC,TFD,CT,CP,MW 
+QD, QR, QT) 


CALL MOISTREM(XW1,XW2,AIR, FLOWRATE, MW, VAIR, DENS, MRR) 
THIS PORTION OF THE PROGRAM PRODUCES THE OUTPUT. 
IF (AIR.EQ.0) THEN 


WRITE (4,12) 
FORMAT (LX, HHH AKKKKKKKKKKEKKKKEKKKKEEEKEEEEKEEEEKRKERKEKEKER 


FHKKKKKKKKKKKKKKKKEKKKKRKRKKKREE! ) 


WRITE (4,15) I 
FORMAT (35X,'CASE ',I2) 


WRITE (4,20) 
FORMAT (LX, TRHERAKKRKKKKKKKEKEKKKKEEKEKKKERKEKERKEKERKERKEREKEEKE 
HFHEKKKKKEKKKKKKEKEKKEKKKKKKKKKEKEKE' ) 


WRITE (4,25) 
FORMAT (25X,'"HYPERBARIC CHAMBER SIMULATION" ') 


WRITE (4,30) 

g 
FORMAT (1X, KKK ERERREKERRRERESE 
FHEKKKKKKEKKEKKEKKEKKEKKEKKEKEKEKEKEKEE ' ) 


WRITE (4,35) 
FORMAT (5X,'INPUT CONDITIONS: ') 


WRITE (4,40) PRESS 
FORMAT (5X, 'PRESSURE (PSIA):',24X,1F6.1) 


WRITE (4,45) TFI,TFD 
FORMAT (5X, 'INITIAL/DESIRED TEMP (DEG F):',9X,1F4.1,'/', 
+1F4.1) 


WRITE (4,50) RHI,RHD 
FORMAT (5X, 'INITIAL/DESIRED REL HUM (%):',10X,1F4.1,'/', 
+1 F 4 ) 


WRITE (4,55) NUMDIV,VAR1 
FORMAT (5X,'# OF DIVERS/WORKING CONDITIONS:',7X,1F4.1, 
+'/',1F2.0) 


WRITE (4,60) CHVOL 
FORMAT (5X, 'CHAMBER VOLUME (CU FT):',17X,1F6.1) 
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WRITE (4,70) FLOWRATE 
FORMAT (5X,'FLOW RATE (CU FT/MIN):',18X,1F6.1) 


WRITE (4,75) DUR 
FORMAT (5X,'MISSION DURATION (HRS):',17X,1F6.1) 


WRITE (4,100) 
FORMAT (1X, VHRR 
tKEKKEKKEKEKKKKEKRKKKRKRKKKRRKKKRKRKRRKEK ) 


WRITE (4,105) 
FORMAT (5X, 'OUTPUT: ') 


WRITE (4,110) PERHE, PERO2 
FORMAT (5X, 'DIVING-GAS MIXTURE (%HE/%02):',9X,1F4.1, 
+'/' 1F4,1) 


WRITE (4,115) DMIN, DMAX 
FORMAT (5X,'DIVING-GAS MIXTURE RANGE (FSW):',7X,1F6.1, 
+'/' 1F6.1) 


WRITE (4,120) O2REQ 
FORMAT (5X, 'OXYGEN REQUIRED (SCFM):',18X,1F8.4) 


WRITE (4,125) LBCO2 
FORMAT (5X,'CO2 PRODUCED (LB):',22X,1F6.1) 


WRITE (4,130) LBS 
FORMAT (5X, 'SODASORB REQUIRED: ',22X,1F6.1) 


WRITE (4,135) TIME 
FORMAT (5X,'TIME TO 95% STEADY STATE (HRS):',9X,1F6.1) 


WRITE (4,140) MW 
FORMAT (5X, 'MOLECULAR WEIGHT (LB/MOLE):',13X,1F8.3) 


WRITE (4,145) DENS 
FORMAT (5X,'GAS MIXTURE DENSITY (LB/CU FT):',9X,1F10.5) 


WRITE (4,150) CT 
FORMAT (5X, 'MAXIMUM COIL TEMPERATURE (DEG F):',7X,1F6.1) 


WRITE (4,155) QD 
FORMAT (5X, 'DEHUMIDIFICATION ENERGY (BTU/MIN):',6xX, 
+1F6.1) 


WRITE (4,160) QR 
FORMAT (5X,'REHEAT ENERGY (BTU/MIN):',16X,1F6.1) 


WRITE (4,165) OT 
FORMAT (5X,'ENERGY REQUIREMENTS (BTU/MIN) :',10X,1F6.1) 
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WRITE (4,170) MRR 
FORMAT (5X, 'MOISTURE REMOVAL RATE (LB/DAY):',9X,1F6.1) 


WRITE (4,175) 
FORMAT (1X, 'X#RRRRRKKKKRKKKKKKKKKKKKKKKKREKRRKEKKEKKKE KKK 


HHEKKKKKKKKKKEKKKKEKKEKKKKKKKEKKRKE ' / 


ELSE 


WRITE (4,200) 
FORMAT CLX, SHREK KKEKKEKEKKKEKEKKEEKKEKEKKEKEKKKKKKEKRKKKKKEKKEKER 
HHEKKKKKKKKKKKKKKKKKKKKKRKKKKEKE!' )/ 


WRITE (4,205) 
FORMAT (25X,'"HYPERBARIC CHAMBER SIMULATION" ') 


WRITE (4,210) 
FORMAT (1X, 1 RR RRKKRRRRKRKRKEKKKKEKKKKKREKRERKEKKKRKKKKKKK KEK 


HHKKKKKKKKKKKKKKKKKKKKKKKKKKK!') 


WRITE (4,215) 
FORMAT (5X, 'INPUT CONDITIONS: ') 


WRITE (4,220) PRESS 
FORMAT (5X,'PRESSURE (PSIA):',24X,1F6.1) 


WRITE (4,225) TFI,TFD 
FORMAT (5X, 'INITIAL/DESIRED TEMP (DEG F):',9X,1F4.1,'/', 
+1F4 J) 


WRITE (4,230) RHI,RHD 
FORMAT (5X, 'INITIAL/DESIRED REL HUM (%):',10X,1F4.1,'/', 
+1F4.1) 


WRITE (4,235) NUMDIV,VAR1 
FORMAT (5X,'# OF DIVERS/WORKING CONDITIONS:',7X,1F4.1, 
+'/' 1F2.0) 


WRITE (4,240) CHVOL 
FORMAT (5X, 'CHAMBER VOLUME (CU FT):',17X,1F6.1) 


WRITE (4,245) 


FORMAT (1X, ' ®RHRK RK KKK KKK KKK KKK KHER KAKA KKK 
FHKKKEKKKKKEKKKKRKKKKKKKKEKKKKKKEEK') 


WRITE (4,250) 
FORMAT (5X, 'OUTPUT: ') 


WRITE (4,255) 
FORMAT (5X, 'DIVING-GAS MIXTURE:',23X,'AIR') 
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WRITE (4,260) 
260 FORMAT (5X,'DIVING-GAS MIXTURE RANGE (FSW):',10X, 
+'0.0/60.0') 


WRITE (4,265) VAIR 
265 FORMAT (5X,'FREE VENTILATION RATE AT DEPTH (CFM):',3X, 
+1F6.1) 


WRITE (4,270) SURAIR 
270 FORMAT (5X,'FREE VENTILATION RATE (SCFM) :',11X,1F6.1) 


WRITE (4,275) TIME 
275 FORMAT (5X,'TIME TO 95% STEADY-STATE (HRS):',9X,1F6.1) 


WRITE (4,280) MW 
280 FORMAT (5X, 'MOLECULAR WEIGHT (LB/MOLE):',15X,1F4.1) 


WRITE (4,285) DENS 
285 FORMAT (5X,'DENSITY OF AIR (LB/CU FT):',14X,1F10.5) 


WRITE (4,290) CT 
290 FORMAT (5X, 'MAXIMUM COIL TEMPERATURE (DEG F):',7X,1F6.1) 


WRITE (4,295) QD 
295 FORMAT (5X,'DEHUMIDIFICATION ENERGY (BTU/MIN):',6X, 
+1F6.1) 


WRITE (4,300) QR 
300 FORMAT (5X,'REHEAT ENERGY (BTU/MIN):',16X,1F6.1) 


WRITE (4,305) QT 
305 FORMAT (5X,'ENERGY REQUIREMENTS (BTU/MIN):',10X,1F6.1) 


WRITE (4,310) MRR 
310 FORMAT (5X,'MOISTURE REMOVAL RATE (LB/DAY):',9X,1F6.1) 


WRITE (4,315) 
q 
315 FORMAT (1X, XR KKKKKKKKKKKKKKKKKK KKK KKKKKKKRKEKKKKKEKKKKKE 
EKKKKKKKEKEKEKKERERKKKEKKEKKKEEKK') 


ENDEF 
500 CONTINUE 


END 
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SUBROUTINE CONVERT (DEPTH, PRESS) 


THIS SUBROUTINE CONVERTS A WORKING DEPTH IN FEET TO A 
PRESSURE IN PSIA. 


REAL DEPTH, PRESS 


A WORKING DEPTH IN FEET. 
A CONVERTED PRESSURE IN PSIA. 


DEPTH 
PRESS 


PRESS=( (DEPTH+33.06633) /33.06633) *14.696 


END 


SUBROUTINE GASMIX(PRESS,O2MIN,O2MAX) 


THIS SUBROUTINE CALCULATES AN ACCEPTABLE DIVING GAS 
MIXTURE BASED UPON A MINIMUM OXYGEN PARTIAL PRESSURE OF 
0.16 ATM AND A MAXIMUM OXYGEN PARTIAL PRESSURE OF 0.6 
ATM. 


REAL PRESS, PATM, O1MIN, O2MIN, O1MAX, O2MAX 


PRESS = THE CHAMBER PRESSURE IN PSIA. 
PATM = THE CHAMBER PRESSURE IN ATM. 


O1MIN = MINIMUM ACCEPTABLE OXYGEN FRACTION IN THE 
DIVING-GAS. 

O2MIN = MINIMUM ACCEPTABLE OXYGEN PERCENTAGE IN THE 
DIVING-GAS. 

O1MAX = MAXIMUM ACCEPTABLE OXYGEN FRACTION IN THE 
DIVING-GAS. 

O2MAX = MAXIMUM ACCEPTABLE OXYGEN PERCENTAGE IN THE 
DIVING-GAS. 


PATM=PRESS/14.696 
O1MIN=0.16/PATM 
O2MIN=01MIN*100.0 
O1MAX=0.6/PATM 
O2MAX=01MAX*100.0 


END 
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SUBROUTINE GASRANGE(PERO2 ,DMIN, DMAX, PSIAMIN, PSIAMAX) 


THIS SUBROUTINE CALCULATES A DEPTH AND PRESSURE RANGE 
FOR A GIVEN DIVING-GAS MIXTURE. 


REAL PERO2, PMIN, DMIN, PMAX, DMAX, PSIAMIN, PSIAMAX 


PERO2 = THE PERCENTAGE OF OXYGEN IN THE DIVING-GAS 
MIXTURE. 


PMIN = THE MINIMUM PRESSURE (ATM) BASED UPON AN OXYGEN 
PARTIAL PRESSURE OF 0.16 ATM. 

DMIN = THE MINIMUM DEPTH (FT) FOR A GIVEN DIVING-GAS. 

PMAX = THE MAXIMUM PRESSURE (ATM) BASED UPON AN OXYGEN 


PARTIAL PRESSURE OF 0.6 ATM. 
DMAX = THE MAXIMUM DEPTH (FT) FOR A GIVEN DIVING-GAS. 
PSIAMIN = THE MINIMUM PRESSURE (PSIA) FOR A GIVEN 
DIVING-GAS. 
PSIAMAX = THE MAXIMUM PRESSURE (PSIA) FOR A GIVEN 
DIVING-GAS. 


PMIN=16.0/PERO2 
DMIN=PMIN*33.06633-33.06633 
PMAX=60.0/PERO2 
DMAX=PMAX*33 .06633-33. 06633 
PSIAMIN=14.696*PMIN 
PSIAMAX=14 .696*PMAX 


END 


SUBROUTINE FREERATE(O2CONS,NUMDIV, TRI, PRESS, VAIR, SURAIR) 


THIS SUBROUTINE CALCULATES THE REQUIRED VENTILATION RATE 
FOR THE CHAMBER WITHOUT CARBON DIOXIDE ABSORBERS. 


REAL O2CONS, NUMDIV, TRI, PRESS, VAIR, SURAIR, CO2PROD, 
+MFLOW 


O2CONS = OXYGEN CONSUMPTION IN CUBIC FEET PER MINUTE. 

NUMDIV = THE NUMBER OF DIVERS IN THE CHAMBER. 

TRI = THE INITIAL CHAMBER TEMPERATURE IN DEGREES 
RANKINE. 

PRESS = THE CHAMBER PRESSURE IN PSIA. 

VAIR = THE VENTILATION RATE AT DEPTH (CU FT/MIN). 

SURAIR = THE VENTILATION RATE AT THE SURFACE (SCFM). 

CO2PROD = THE TOTAL AMOUNT OF CARBON DIOXIDE PRODUCED 

(CU FT/MIN). 
MFLOW = THE MASS FLOW RATE OF CARBON DIOXIDE (LB/MIN). 
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CO2PROD=0.9*02CONS*NUMDIV 
MFLOW=0. 11372474*CO2PROD 
VAIR=3 .31817426*MFLOW*TRI 
SURAIR=(VAIR*PRESS*529.67) /(14.696*TRI) 


END 


SUBROUTINE ABSRATE(NUMDIV,O2REQ,02CONS, DUR, LBCO2, LBS) 


THIS SUBROUTINE CALCULATES THE REQUIRED VENTILATION RATE 
OF THE CHAMBER IF CARBON DIOXIDE ABSORBERS ARE UTILIZED. 


REAL NUMDIV, O2REQ, O2CONS, DUR, LBCO2, MFLOW, LBS, 
+CO2PROD 


NUMDIV = THE NUMBER OF DIVERS IN THE CHAMBER. 

O2REQ = THE AMOUNT OF PURE OXYGEN REQUIRED (SCFM). 

O2CONS = THE OXYGEN CONSUMPTION (SCFM). 

DUR = THE MISSION DURATION IN HOURS. 

LBCO2 = THE AMOUNT OF CO2 PRODUCED (LB). 

MFLOW = THE MASS FLOWRATE OF CARBON DIOXIDE (LB/MIN). 

CO2PROD = CARBON DIOXIDE PRODUCTION (CU FT/MIN). 

LBS = AMOUNT OF SODASORB REQUIRED FOR CO2 ABSORPTION 
(LBS) . 


O2REQ=NUMDIV*0O2CONS 
CO2PROD=0.9*02CONS*NUMDIV 
MFLOW=0.11372474*CO2PROD 
LBCO2=MFLOW*60.0*DUR 
LBS=LBCO2/0.245 


END 


SUBROUTINE STSTATE(CHVOL, O02CONS,NUMDIV, TRI, TIME) 


THIS SUBROUTINE CALCULATES THE TIME REQUIRED FOR THE 
CHAMBER TO REACH A 95% STEADY-STATE CONDITION. 


REAL CHVOL, O2CONS, NUMDIV, TRI, TIME, CO2PROD, MFLOW, 
+VAIR 


CHVOL = THE INTERIOR CHAMBER VOLUME (CUBIC FEET). 
O2CONS = OXYGEN CONSUMPTION (CU FT/MIN). 
NUMDIV = THE NUMBER OF DIVERS IN THE CHAMBER. 
TRI = THE INITIAL TEMPERATURE OF THE CHAMBER IN DEGREES 
RANKINE. 
TIME = THE TIME REQUIRED TO REACH 99% STEADY-STATE 
(HRS) . 
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CO2PROD = THE CARBON DIOXIDE PRODUCTION (CU FT/MIN). 
MFLOW = THE MASS FLOWRATE OF CARBON DIOXIDE (LB/MIN). 
VAIR = THE VENTILATION RATE AT DEPTH (CU FT/MIN). 


CO2PROD=0.9*02CONS*NUMDIV 

MFLOW=0.11372474*CO2PROD 

VAIR=3.31817426*MFLOW*TRI 

TIME=(-0.0167*CHVOL/VAIR) *LOG(1.0-(0.2863*VAIR) 
+  /(MFLOW*TRI) ) 


END 
SUBROUTINE MOLWT(PERO2,MW,AIR) 


THIS SUBROUTINE CALCULATES THE MOLECULAR WEIGHT OF THE 
DIVING-GAS MIXTURE. 


REAL PERO2, MW, X02, XHE, AIR 
PERO2 = THE PERCENTAGE OF OXYGEN IN THE DIVING-GAS 


MIXTURE. 
MW = THE MOLECULAR WEIGHT OF THE MIXTURE (LB/LB-MOLE). 


XO2 = THE MOLE FRACTION OF OXYGEN IN THE MIXTURE. 
XHE = THE MOLE FRACTION OF HELIUM IN THE MIXTURE. 
AIR = A VARIABLE TO DETERMINE IF THE DIVING-GAS IS AIR. 


IF (AIR.EQ.1) THEN 
MW=29.0 
ELSE 
XO2=PERO2/100.0 
XHE=1.0-X02 
MW= (X02*32.00)+(XHE*4.0026) 
ENDIF 
END 


SUBROUTINE SPECHEAT (PERO2,CP,AIR) 


THIS SUBROUTINE CALCULATES THE SPECIFIC HEAT OF THE 
DIVING-GAS MIXTURE. 


REAL PERO2, CP, X02, XHE, AIR 


PERO2 = THE PERCENTAGE OF OXYGEN IN THE DIVING-GAS 
MIXTURE. 


105 





Le eke ae: 


Orme) O OWA OO] a ae 


CP = THE SPECIFIC HEAT OF THE DIVING-GAS MIXTURE 
(BTU/MOLE) . 


XO2 = THE MOLE FRACTION OF OXYGEN IN THE MIXTURE. 
XHE = THE MOLE FRACTION OF HELIUM IN THE MIXTURE. 
AIR = A VARIABLE TO DETERMINE IF THE DIVING-GAS IS AIR. 


IF (AIR.EQ.1) THEN 
CP=6.95 
ELSE 


XO2=PERO2/100.0 
XHE=1.0-X02 
CP=(X02*6.95)+(XHE*4.965) 


ENDIF 


END 


SUBROUTINE DENSITY (PERO2,TRI,PRESS,MW,DENS, BC, BF,AIR) 


THIS SUBROUTINE CALCULATES THE DENSITY OF THE DIVING-GAS 
MIXTURE. 


REAL PERO2, TRI, PRESS, MW, DENS, BC, BF, CONST, MOLVOL, 
+AIR 


PERO2 = THE PERCENTAGE OF OXYGEN IN THE DIVING-GAS 
MIXTURE. 
TRI = THE INITIAL CHAMBER TEMPERATURE IN DEGREES 
RANKINE. 
PRESS= THE INITIAL CHAMBER PRESSURE IN PSIA. 
MW = THE MOLECULAR WEIGHT OF THE DIVING-GAS MIXTURE 
(LB/LB-MOLE) . 
DENS = THE DENSITY OF THE DIVING-GAS MIXTURE (LB/CU 
FEET) . 
BC = THE SECOND VIRIAL COEFFICIENT (CU CM/MOLE). 
BF = THE SECOND VIRIAL COEFFICIENT (CU FT/MOLE) . 
CONST = RT/P 
MOLVOL = MOLAR VOLUME (CU FT/MOLE) 
AIR = A VARIABLE TO DETERMINE IF THE DIVING-GAS IS AIR. 


IF (AIR.EQ.1) THEN 
DENS= (MW*PRESS*144.0) /(1544.0*TRI) 
ELSE 


XO2=PERO2/100.0 
XHE=1.0-X02 
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BC=(11.44*XHE**2) -(8.8*XHE*XO2) -(16.24*X02**2) 
BF=BC*3 .5315E-05*454.0 
CONST=(1544.0*TRI) / (PRESS*144.0) 

MOLVOL= (CONST+SQRT( (CONST**2)+(4.0*CONST*BF) ))/2.0 
DENS=MW /MOLVOL 


ENDIF 


END 


SUBROUTINE PSYCHART(TFI,TRI, PRESS, XWS,H) 


THIS SUBROUTINE CALCULATES THE MOISTURE MOLE FRACTION 
AND ENTHALPY OF MOISTURE UTILIZING THE PSYCHROMETRIC 
CHART. 


REAL TFI, TRI, PRESS, XWS, H, PS, DPS, VL 


TFI = THE INITIAL CHAMBER TEMPERATURE IN DEGREES 

FARENHEIT. 

THE INITIAL CHAMBER TEMPERATURE IN DEGREES 

RANKINE. 

PRESS = THE INITIAL CHAMBER PRESSURE IN PSIA. 

XWS = THE MOISTURE MOLE FRACTION (MOLES WATER VAPOR/MOLE 
DRY GAS). 


TRI 


H = THE ENTHALPY OF MOISTURE CONTENT (BTU/MOLE DRY GAS). 
PS = THE SATURATION PRESSURE OF WATER (PSIA). 
DPS = THE CHANGE IN SATURATION PRESSURE OF WATER FOR 


PRESSURES ABOVE 100 PSIA. 
VL = THE SPECIFIC VOLUME OF THE SATURATED LIQUID (CU 
FT/LB) . 


THIS PORTION OF THE SUBROUTINE CALCULATES THE SATURATION 
PRESSURE OF WATER GIVEN THE DRY BULB TEMPERATURE. 


PS=(1.9285E-11*TFI**5)+(9.8683E-10*TFI**4) + 
+(3.3193E-07*TFI**3)+(1.9916E-05*TFI**2) + 
+(1.0488E-03*TFI) +2.2080E-02 


THIS PORTION OF THE SUBROUTINE CORRECTS THE SATURATION 
PRESSURE FOR PRESSURES ABOVE 100 PSIA USING POYNTING'S 
EQUATION. 


IF (PRESS.GT.100.0) THEN 
IF (TFI.GT.32.0.AND.TFI.LE.52.0) THEN 
VL=0.01602 
ELSEIF (TFI.GT.52.0.AND.TFI.LE.60.0) THEN 
VL=0.01603 
ELSEIF (TFI.GT.60.0.AND.TFI.LE.66.0) THEN 
VL=0.01604 
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ELSEIF (TFI.GT.66.0.AND.TFI.LE.71.0) THEN 


VL=0.01605 

ELSEIF (TFI.GT.71.0.AND.TFI.LE.76.0) THEN 
VL=0. 01606 

ELSEIF (TFI.GT.76.0.AND.TFI.LE.80.0) THEN 
VL=0.01607 

ELSEIF (TFI.GT.80.0.AND.TFI.LE.84.0) THEN 
VL=0.01608 

ELSEIF (TFI.GT.84.0.AND.TFI.LE.87.0) THEN 
VL=0.01609 

ELSEIF (TFI.GT.87.0.AND.TFI.LE.91.0) THEN 
VL=0.01610 

ELSEIF (TFI.GT.91.0.AND.TFI.LE.94.0) THEN 
VL=0.01611 

ELSEIF (TFI.GT.94.0.AND.TFI.LE.97.0) THEN 
VL=0.01612 

ELSEIF (TFI.GT.97.0.AND.TFI.LE.100.0) THEN 
VL=0.01613 

ELSEIF (TFI.GT.100.0) THEN 
VL=0.01614 

ENDIF 


DPS=(1.67876*PS*VL*PRESS) /TRI 
PS=PS+DPS 


GO TO 40 
ELSE 
GO TO 40 
ENDIF 


XWS=PS/ (PRESS-PS) 
H=20000.0*XWS 


END 


SUBROUTINE XWSH(RHI,RHD, XWS, XW1, XW2,H,H1,H2) 


THIS SUBROUTINE CALCULATES THE INITIAL AND DESIRED 
MOISTURE MOLE FRACTIONS AND MOLAR ENTHALPIES. 


REAL RHI, RHD, XWS, XW1, XW2, H, H1, H2, RH1, RH2 


RHI = THE INITIAL RELATIVE HUMIDITY IN THE CHAMBER IN 


PERCENT. 

RHD = THE DESIRED RELATIVE HUMIDITY IN THE CHAMBER IN 
PERCENT. 

XWS = THE MOISTURE MOLE FRACTION FOR WATER VAPOR AT 
SATURATION. 
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AW1 = THE MOISTURE MOLE FRACTION AT THE INITIAL RELATIVE 
HUMIDITY. 
XW2 = THE MOISTURE MOLE FRACTION AT THE DESIRED RELATIVE 
HUMIDITY. 
H = THE MOLAR ENTHALPY OF MOISTURE CONTENT AT 
SATURATION. 


Hl = THE MOLAR ENTHALPY OF MOISTURE CONTENT AT THE 
INITIAL R.H. 

H2 = THE MOLAR ENTHALPY OF MOISTURE CONTENT AN THE 
DESIRED R.H. 

RH1 = THE INITIAL RELATIVE HUMIDITY IN THE CHAMBER 
(FRACTION) . 

RH2 = THE DESIRED RELATIVE HUMIDITY IN THE CHAMBER 
(FRACTION) . 


RH1=RHI/100.0 
RH2=RHD/100.0 
XW1=RH1*XWS 
XW2=RH2*XWS 
H1=RH1*H 
H2=RH2*H 


END 
SUBROUTINE COILTEMP (XW2, PRESS, CT) 


THIS SUBROUTINE CALCULATES THE COIL TEMPERATURE REQUIRED 
FOR DEHUMIDIFICATION TO A DESIRED RELATIVE HUMIDITY. 


REAL XW2, PRESS, PS, CT 
XW2 = THE MOISTURE MOLE FRACTION AT THE DESIRED RELATIVE 
HUMIDITY. 

PRESS = THE INITIAL CHAMBER PRESSURE (PSIA). 

PS = THE SATURATION PRESSURE OF WATER (PSIA). 

CT = THE COIL TEMPERATURE REQUIRED FOR DEHUMIDIFICATION. 
PS=(XW2*PRESS) /(1.0+XW2) 

CT=(2.4135E02*PS**5) -(8.4538E02*PS**4) + 
+(1.1635E03*PS**3) -(8.1910E02*PS**2) + 

+(3.5407E02*PS) +7.0106 


END 


109 





AQAA 


YO AWIY A OOOO Ga Ga a 


QAANAANA 


7 OQ) Ae 


SUBROUTINE MGMENTH(TFI,CT,CP,H1,H2,MW,HA,HC) 
THIS SUBROUTINE CALCULATES THE MOIST-GAS-MIXTURE 
ENTHALPIES AT THE INITIAL CONDITIONS AND COIL 
TEMPERATURE. 

REAL TFI, CT, CP, H1, H2, MW, HA, HC 


TFI = THE INITIAL CHAMBER TEMPERATURE IN DEGREES 


FARENHEIT. 

CT = THE REQUIRED COIL TEMPERATURE FOR DEHUMIDIFICATION 
(F). 

CP = THE SPECIFIC HEAT OF THE GAS MIXTURE (BTU/MOLE) . 

H1 = THE MOLAR ENTHALPY OF MOISTURE CONTENT AT THE 
INITIAL R.H. 

H2 = THE MOLAR ENTHALPY OF MOISTURE CONTENT AT THE 
DESIRED R.H. 

MW = THE MOLECULAR WEIGHT OF THE MIXTURE (LB/LB-MOLE) . 

HA = THE MOIST-GAS-MIXTURE ENTHALPY AT THE INITIAL 
CONDITIONS. 

HC = THE MOIST-GAS-MIXTURE ENTHALPY AT THE DESIRED 
CONDITIONS. 


HA=((CP*TFI) +H1) /MW 
HC=((CP*CT) +H2) /MW 


END 


SUBROUTINE ENERGY (DENS, AIR, VAIR, FLOWRATE,HA,HC,TFD,CT, 
+CP,MW,QD,QR, QT) 


THIS SUBROUTINE CALCULATES THE ENERGY REQUIRED FOR 
DEHUMIDIFICATION AND REHEAT. 


REAL DENS, AIR, VAIR, FLOWRATE, HA, HC, TFD, CT, CP, MW, 
+OD, DH OR, OT 


DENS = THE DENSITY OF THE DIVING-GAS MIXTURE (LB/CU FT). 

AIR = A VARIABLE TO DETERMINE IF THE DIVING-GAS IS AIR. 

VAIR = THE VENTILATION OR CIRCULATION RATE AT DEPTH (CU 
FT/MIN). 

FLOWRATE = THE CIRCULATION RATE (CU FT/MIN). 


HA = THE MOIST-GAS-MIXTURE ENTHALPY AT THE INITIAL 
CONDITIONS. 

HC = THE MOIST-GAS-MIXTURE ENTHALPY AT THE DESIRED 
CONDITIONS. 


TFD = THE DESIRED CHAMBER TEMPERATURE AFTER REHEAT (F). 
CT = THE COIL TEMPERATURE REQUIRED FOR DEHUMIDIFICATION 
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CP = THE SPECIFIC HEAT OF THE DIVING-GAS MIXTURE 
(BTU/MOLE) . 
MW = THE MOLECULAR WEIGHT OF THE MIXTURE (LB/LB-MOLE) . 


QD = THE ENERGY REQUIRED FOR DEHUMIDIFICATION (BTU/MIN) . 

DH = THE CHANGE IN ENTHALPY FROM THE COIL TEMPERATURE TO 
THE DESIRED CHAMBER TEMPERATURE (BTU/MIN) . 

QR = THE ENERGY REQUIRED FOR REHEAT (BTU/MIN). 

QT = THE TOTAL ENERGY REQUIRED DEHUMIDIFICATION AND 


REHEAT (BTU/MIN) . 


IF (AIR.EQ.0) THEN 
VAIR=FLOWRATE 


ELSE 
VAIR=VAIR 


ENDIF 


QD=DENS*VAIR* (HA-HC) 
DH=CP* (TFD-CT) /MW 
QR=DENS*VAIR*DH 
QT=QD+QR 


END 


SUBROUTINE MOISTREM(XW1,XW2,AIR, FLOWRATE, MW,VAIR, 
+DENS, MRR) 


THIS SUBROUTINE CALCULATES THE MOISTURE REMOVAL RATE 
REQUIRED FOR DEHUMIDIFICATION. 


REAL W1, W2, XW1, XW2, AIR, FLOWRATE, MW, VAIR, DENS, 
+MRR 


Wl = THE HUMIDITY RATIO AT THE INITIAL RELATIVE 
HUMIDITY. 
W2 = THE HUMIDITY RATIO AT THE DESIRED RELATIVE 
HUMIDITY. 
XW1 = THE MOISTURE MOLE FRACTION AT THE INITIAL RELATIVE 
HUMIDITY. 
XW2 = THE MOISTURE MOLE FRACTION AT THE DESIRED RELATIVE 
HUMIDITY. 
AIR = A VARIABLE TO DETERMINE IF THE DIVING-GAS IS AIR. 
FLOWRATE = THE CIRCULATION RATE (CU FT/MIN). 
MW = THE MOLECULAR WEIGHT OF THE DIVING-GAS MIXTURE 
(LB/LB-MOLE) . 
VAIR = THE VENTILATION RATE AT DEPTH (CU FT/MIN). 
DENS = THE DENSITY OF THE DIVING-GAS MIXTURE (LB/CU FT). 
MRR = THE REQUIRED MOISTURE REMOVAL RATE FOR THE 
DEHUMIDIFICATION PROCESS (LB/DAY). 


ieee 





IF (AIR.EQ.0) THEN 
VAIR=FLOWRATE 


ELSE 
VAIR=VAIR 
ENDIF 
W1=XW1* (18.0158/MW) 
W2=XW2*(18.0158/MW) 
MRR= (W1-W2) *VAIR*DENS*1440.0 


END 
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